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METHODS AND SYSTEMS FOR 
CHARACTERIZATION OF AN ANOMALY 
USING INFRARED FLASH 
THERMOGRAPHY 

5 

I. CROSS REFERENCE TO RELATED 
APPLICATIONS 

This application claims the benefit of the filing date of U.S. 
provisional patent application No. 61/258,266, incorporated 
herein by reference, which was filed on Nov. 5, 2009, by the 
same inventor of this application. 

II. ORIGIN OF THE INVENTION 

The invention described herein was made by an employee 
of the United States Government and may be manufactured 
and used by or for the Government of the United States of 
America for governmental purposes without the payment of 
any royalties thereon or therefor. 

III. FIELD OF THE INVENTION 

The present invention generally relates to the field of non- 
destructive evaluation (NDE) and of processing of data 
acquired from a system performing NDE. More particularly, 
the invention relates to methods and systems for analyzing 
and processing digital data comprising a plurality of infrared 
(IR) video images acquired by a system used for non-destruc- 
tive evaluation. 

IV. BACKGROUND OF THE INVENTION 

The embodiments described herein relate generally to the 
field of non-destructive evaluation (NDE). Infrared flash (or 
pulsed) thermography is an example of a technique for NDE 
used primarily in the inspection of thin nomnetallic materials, 
such as laminated or bonded composites in the aerospace 
industry. IR flash thermography is primarily used to detect 
delamination-like anomalies, although surface cracks may 
also be detected to some extent. In one embodiment, a single- 
sided or reflection technique is used wherein the flash lamp 
(heat source) and the IR camera (detector) are on the same 
side of the test object undergoing inspection. 

The hardware equipment for an IR flash thermography 
system comprises a flash lamp (source of light/heat), a flash 
hood, a flash power supply/trigger unit, a flash duration con- 
troller, an IR camera for capturing video images, data acqui- 
sition electronics, and a computer. The computer is used for 
controlling the flash trigger, for acquiring video data from the 
IR camera, for displaying data, and for post-processing of the 
acquired data. 

In one example of an NDE technique using IR flash ther- 
mography, a plate is provided as a test object with a round 
delamination in the center. After applying heat to the top 
surface of the test object by triggering the flash lamp, the top 
surface area surrounding the anomaly cools faster than the top 
surface (footprint) area above the anomaly. The IR camera 
captures a sequence of images of the surface temperature in 
terms of pixel intensity and shows the anomaly as a hot spot 
(e.g., an area warmer than the surrounding area or the refer- 
ence region of interest (ROI)). The hot spot is about the size 
and shape of the anomaly footprint. Relative pixel intensity, 
i.e., the difference in pixel intensity between the hot spot 
(measurement ROI) and surrounding area (reference ROI), 
varies with the post-flash time. Deeper anomalies appear at 
later times in the IR video data compared to the near surface 
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anomalies. After the appearance of an anomaly in the IR 
video data, the relative pixel intensity continues to increase 
with time. The relative pixel intensity of the anomaly reaches 
a peak at a certain time and then the relative pixel intensity 
decays until the temperature of the indication area and the 
temperature of the surrounding area become equal. 

V. SUMMARY 

The exemplary embodiments described herein extend the 
theory and application of IR flash thermography to not just 
anomaly detection, but to anomaly characterization. The 
embodiments described herein examine the relative pixel 
intensity and the relative temperature between the measure- 
ment region of interest (ROI) and the reference ROI by using 
a combination of equations generically defined as the “nor- 
malized pixel intensity contrast” and the “normalized tem- 
perature contrast.” By calculating a “normalized pixel inten- 
sity contrast” during the processing of the acquired IR video 
data, the embodiments comprise the development of a simu- 
lation model, a data analysis approach and method, and cor- 
responding software for implementing the model, approach, 
and method with a computer, the methods being collectively 
referred to as the “IR Contrast” or the “Koshti Contrast” 
method. The software for implementing the “IR Contrast” 
method on a computer comprises an IR contrast prediction 
model that employs contrast matching algorithms for post- 
processing of the data acquired by an IR camera operating as 
part of a flash thermography system or as part of other non- 
destructive evaluation systems. 

The IR Contrast method is focused on anomaly character- 
ization as opposed to mere anomaly detection. The approach 
for the IR Contrast method relates thermal measurements, 
such as the contrast evolution (or contrast signal) amplitude 
of an anomaly, in terms of the normalized anomaly contrast as 
a function of the frame number or the post-flash time. Under 
some constraints, the normalized anomaly contrast is mainly 
influenced by factors “related” to an anomaly such as the 
material thermal properties (diffusivity and emissivity), 
anomaly size, the depth, and the gap thickness (heat trans- 
missivity or thermal mismatch factor). The contrast is also 
affected by factors that are “unrelated” to the anomaly such as 
the part curvature, uneven part cross section, varying thermal 
properties along the part surface, surface cracks, surface tex- 
ture, uneven flash heating, flash afterglow, camera pixel size, 
focusing and proximity to the part edge or the flash boundary. 
The IR Contrast method seeks to either minimize the effects 
of the unrelated factors in order to improve accuracy of the IR 
Contrast data analysis or to control the unrelated factors to 
improve repeatability of the IR Contrast data analysis. Sev- 
eral amplitude and time-related thermal measurements (e.g. 
peak contrast and peak contrast time, peak product time, 
persistence energy time, etc.) are computed from the mea- 
sured contrast evolution. These parameters bear correlation to 
the anomaly depth, width, and gap thickness (or heat trans- 
missivity through the anomaly). In other embodiments 
described herein, the approach of the IR Contrast method 
comprises measurement of the width of the indication using 
the half-max technique, as will be explained in further detail 
later. 

The IR Contrast method also comprises steps for the cali- 
brated simulation of the contrast evolution. The software 
calculates the diameter and depth as inputs, called equivalent 
flat bottom hole (EFBH) inputs, to the simulation that provide 
the best match to the measured contrast evolution. Thus, the 
anomaly depth and width are mapped into an EFBH diameter 
and depth. If the ratio of an EFBH diameter to half-max 
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diameter is close to one for a round anomaly, then the 
anomaly is similarto a FBH or a gapping round delamination. 
Similarly, if the ratio of an equivalent uniform gap (EUG) 
diameter to half-max diameter is close to one for a round 
anomaly, then the anomaly is similar to a round void with a 5 
uniform gap thickness. The IR Contrast method hypothesizes 
that the EFBH depth is an estimation of the anomaly depth if 
the anomaly is similar to a gapping delamination or a FBF1. 
Similarly, the EUG depth is an estimation of the anomaly 
depth if the anomaly is similar to a void with a uniform gap 10 
thickness. Embodiments of the IR Contrast method allow 
characterization of an anomaly in tenns of its depth, width, 
and gap thickness. In another embodiment, the software 
implementing the IR Contrast method is programmed for 15 
variously different materials and therefore provides a quick 
approximate prediction of the contrast on a given delaminat- 
ing-like anomaly. The software can also be used as a training 
tool for the application of flash thermography or other equiva- 
lent non-destructive evaluation systems. 20 

The computer-implemented IR Contrast software com- 
prises two methodologies to analyze the flash thermography 
contrast data for characterization of delaminating-like 
anomalies . The first methodology uses measurements such as 
the peak contrast, persistence energy time, withholding time, 25 
and peak time, which can be plotted against the location of the 
measurement. These measurements show some correlation to 
the anomaly depth, width, and gap thickness. The second 
methodology comprises matching the calibration-based pre- 
dicted contrast evolution to the measured contrast evolution. ,( 1 
The matching of the simulated contrast evolution to the mea- 
sured contrast provides a size and depth based on the equiva- 
lent flat bottom hole (EFBH) or equivalent uniform gap 
(EUG) with known attenuation. In other words, the EFBH (or , 
EUG) depth and diameter are correlated to the width and 
depth of the measured anomaly. 

Another method known as "half-max edge detection” 
complements the IR Contrast EFBH (or EUG) width predic- 
tions, and the reverse is also true. The difference in the two 40 
width measurements is considered to be due to non-unifor- 
mity of the gap thickness, value of gap thickness, and length 
effect as well as to the limitations of the two estimations. The 
embodiments described herein provide information on the IR 
Contrast/half-max analysis of the flash thermography data. 45 
Comparison of the adjusted EUG diameter and EFBH diam- 
eter with the half-max width allows assessment of the gap 
thickness. Comparison of the peak product time with the 
square of the half-max width also allows assessment of the 
gap thickness. 50 

Other aspects and advantages of the embodiments 
described herein will become apparent from the following 
description and the accompanying drawings, illustrating the 
principles of the embodiments by way of example only. 

55 

VI. BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a schematic diagram of an exemplary single-sided 
flash thermography system in accordance with embodiments 
described herein. 60 

FIG. 2(a) is a schematic diagram of a plate depicting a 
gapping delamination with measurement and reference 
regions of interest (ROIs) in accordance with embodiments 
described herein. 

FIG. 2(b) is a single infrared (IR) video image of a flat 65 
bottom hole in a test object in accordance with some embodi- 
ments. 
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FIG. 3 is a graph of the normalized temperature contrast as 
a function of the normalized time in accordance with embodi- 
ments described herein. 

FIG. 4 is a graph of the normalized reference contrast 
evolution curve in accordance with embodiments described 
herein. 

FIG. 5 is a graph of the flaw size parameter as a function of 
the depth and diameter in accordance with embodiments 
described herein. 

FIG. 6 is a graph of the flaw size parameter as a function of 
gamma in accordance with embodiments described herein. 

FIG. 7 is a graph of the transformation coefficient CTas a 
function of the flaw size parameter in accordance with 
embodiments described herein. 

FIG. 8 is a graph of the inverse of the transformation 
coefficient as a function of gamma in accordance with 
embodiments described herein. 

FIG. 9(a) is a graph of the diameter constant and amplitude 
constant calibration curves. FIG. 9(b) shows an example of 
calibration curves for no attenuation (Flat Bottom Hole) and 
with 30 percent and 50 percent attenuation in accordance with 
embodiments described herein. 

FIG. 10 is a flow diagram illustrating the process for simu- 
lating the contrast evolution. 

FIG. 11 is a flow diagram illustrating the direct estimation 
of the equivalent flat bottom hole depth and diameter in accor- 
dance with embodiments described herein. 

FIG. 12 is a table showing an example of calculation inputs 
in direct depth estimation. 

FIG. 13 is a graph of the measured raw contrast evolution, 
the filtered contrast evolution, and the simulated (predicted) 
contrast evolution in accordance with embodiments 
described herein. 

FIG. 14 is a flow diagram showing the process of contrast 
simulation match with the measured contrast evolution in 
accordance with embodiments described herein. 

FIG. 15 is a frame from a raw data-cube with temperature 
versus time curves at selected regions of interest (ROI) 
including a reference region of interest. 

FIG. 16 is a frame from the contrast data-cube with the 
contrast versus time curves at selected region of interests in 
accordance with embodiments described herein. 

FIG. 17 is a flow diagram of how the measured flash ther- 
mography data is processed with four processing outcomes 
including: (1 ) the normalized datacube and feature measure- 
ment (e.g. peak contrast, peak time) on the entire contrast data 
cube and plotting of feature frame images similar to pixel 
intensity frame images; (2) the half-max anomaly bound- 
aries; (3) measurements on the contrast evolution that relate 
to flaw depth and gap thickness; and (4) simulation match 
equivalent depth and diameter measurement and width ratio 
estimation that relates to the gap thickness, in accordance 
with embodiments described herein. 

FIG. 18 is a graph of the measurement and reference 
regions of interest decay curves plotted in log time in accor- 
dance with embodiments described herein. 

FIG. 19 is a graph showing the peak contrast and the peak 
contrast time in accordance with embodiments described 
herein. 

FIG. 20 is a graph of the peak contrast and peak contrast 
time simulation (contrast map) using the ThermoCalc simu- 
lation. 

FIG. 21(a) and FIG. 21 (b) are graphs of the peak contrast 
and peak time simulation (contrast map) using IR Contrast 
calibration on a reinforced carbon-carbon part in accordance 
with embodiments described herein. 
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FIG. 22 illustrates the use of the adjusted flaw size param- 
eter in prediction of the peak contrast in accordance with 
embodiments described herein. 

FIG. 23(a) shows the back side of a test plate with flat 
bottom holes and slots as well as a cropped area from a 5 
selected frame from the flash thermography data. Two iden- 
tical IR images that include the region of holes numbered 3, 4, 
and 5 are provided; one without traces and one with half-max 
boundary traces, measurement ROI trace, and reference ROI 
trace, in accordance with embodiments described herein, to 
FIG. 23(6) shows both the actual width and the equivalent 
diameter on holes numbered 3, 4, and 5 in a reinforced car- 
bon-carbon (RCC) test specimen in accordance with embodi- 
ments described herein. 

FIG. 24 shows an early IR frame of the flash thermography 1 5 
data for a graphite epoxy laminate. The image shows the 
weave like pattern called texture and the simple contrast 
evolutions illustrating the positive and negative texture peaks 
in accordance with embodiments described herein. 

FIG. 25 displays multiple graphs illustrating the effect of 20 
the frame rate on the contrast evolution of shallow anomalies 
in accordance with embodiments described herein. 

FIG. 26 is a flow diagram of the process of measuring the 
thermal diffusivity of the test piece with flat bottom holes in 
accordance with embodiments described herein. 25 

FIG. 27 displays multiple graphs illustrating the effect of 
the diffusivity and depth constant on the contrast evolution in 
accordance with embodiments described herein. 

FIG. 28 provides multiple graphs illustrating the effect of 
change in the diameter on the contrast evolution in accor- 30 
dance with embodiments described herein. 

FIG. 29 illustrates the effect of change in the attenuation p 
on the contrast evolution in accordance with embodiments 
described herein. 

FIG. 30 is a flow diagram for generating calibration curves 35 
for flat bottom holes (no attenuation), calculating the adjusted 
flaw size parameter as a function of gamma, and creating a 
contrast map in accordance with embodiments described 
herein. 

FIG. 31(a) shows results of a ThermoCalc simulation to 40 
generate relationships between the peak time ratio and the 
peak contrast ratio for various gamma values. The peak time 
ratio for the 50% attenuation is read from FIG. 31(a) and the 
results are shown in FIG. 31(6). 

FIG. 32 is a flow diagram for the process of generating 45 
calibration curves for uniform gap delamination-like voids 
with selected attenuation and of creating contrast maps in 
accordance with embodiments described herein. 

FIG. 33(a) shows experimental correlation between the 
width ratio and gap thickness in RCC material. FIG. 33(6) 50 
shows experimental correlation between the gap correlation 
parameter defined as peak product time divided by the square 
of the half-max width and the gap thickness in RCC material 
in accordance with embodiments described herein. 

FIG. 34 shows a ThermoCalc simulated data correlation 55 
between the gap correlation parameter defined as peak prod- 
uct time divided by 2.3 power of the simulated void width and 
the simulated void gap thickness in a material. 

FIG. 35 is a block diagram illustrating an apparatus for 
characterizing anomalies in a material in accordance with 60 
embodiments described herein. 

While the embodiments described herein are subject to 
various modifications and alternative forms, the specific 
embodiments are shown by way of example in the drawings 
and the accompanying detailed description. It should be 65 
understood, however, that the drawings and detailed descrip- 
tion are not intended to limit the invention to the particular 
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embodiments described herein. This description is instead 
intended to cover all modifications and alternatives as falling 
within the scope of the invention as defined by the appended 
claims. 

VII. DETAILED DESCRIPTION 

The objects in the drawings are not necessarily to scale and 
certain features may be shown exaggerated in scale or in 
somewhat generalized or schematic form in the interest of 
clarity and conciseness. In the description which follows like 
parts may be marked throughout the specification and draw- 
ings with the same reference numerals. The foregoing brief 
description of the drawings is provided for a more complete 
understanding thereof. It should be understood, however, that 
the embodiments described are not limited to the precise 
arrangements and configurations shown. Although the design 
and use of various embodiments are discussed in detail below, 
it should be appreciated that the embodiments described rep- 
resent concepts that may be embodied in a wide variety of 
contexts. The specific aspects and embodiments described 
herein are merely illustrative, for it would be impossible or 
impractical to include all of the possible embodiments and 
contexts in this detailed description. Upon reading this 
description, alternative embodiments will be apparent to per- 
sons of ordinary skill in the art. 

The embodiments described herein examine the relative 
pixel intensity and the relative temperature between the mea- 
surement region of interest (ROI) and the reference ROI by 
using a combination of equations generically defined as the 
“normalized pixel intensity contrast” and the “normalized 
temperature contrast.” By calculating a “normalized pixel 
intensity contrast” during the processing of the IR video data 
acquired, the embodiments comprise the development of a 
simulation model, a data analysis approach and method, and 
corresponding software for implementing the model, 
approach, and method with a computer, collectively referred 
to as the “IR Contrast” or the “Koshti Contrast” method. The 
software for implementing the “IR Contrast” method on a 
computer comprises an IR contrast prediction model that 
employs contrast matching algorithms for post-processing of 
the data acquired from an IR camera, with the camera oper- 
ating in some embodiments as part of a flash thermography 
system or in other embodiments as part of other non-destruc- 
tive evaluation systems. 

The IR Contrast method is geared toward anomaly charac- 
terization as opposed to mere anomaly detection. The 
approach for the IR Contrast method relates thermal measure- 
ments, such as the contrast evolution (or contrast signal) 
amplitude of an anomaly, in terms of the normalized anomaly 
contrast as a function of the frame number or the post-flash 
time. Under some constraints, the normalized anomaly con- 
trast is mainly influenced by factors “related” to an anomaly 
such as the material thermal properties (diffusivity and emis- 
sivity), the anomaly size, the anomaly depth, the gap thick- 
ness, and the heat transmissivity or thermal mismatch factor. 
The contrast is also affected by factors that are “unrelated” to 
the anomaly such as the part curvature, uneven part cross 
section, varying thermal properties along the part surface, 
surface cracks, surface texture, uneven flash heating, flash 
afterglow, camera pixel size, focusing, and proximity to the 
part edge or the flash boundary. The IR Contrast method seeks 
to either minimize the effect of the unrelated factors in order 
to improve accuracy of the IR Contrast analysis or to control 
the unrelated factors to improve repeatability of the IR Con- 
trast analysis. Several amplitude and time-related thermal 
measurements (e.g. peak contrast and peak contrast time, 
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peak product time, persistence energy time, etc.) are com- 
puted from the measured contrast evolution. These param- 
eters bear some correlation to the anomaly depth, width, and 
gap thickness (or heat transmissivity through the anomaly). 
The approach also involves measurement of the indication 5 
using the half-max technique. 

The IR Contrast method also comprises steps for calibrated 
simulation of the contrast evolution. The software calculates 
the diameter and depth as input iteratively, called equivalent 
flat bottom hole (EFBH) inputs, to the simulation that pro- to 
vides the best match to the measured contrast evolution. Thus, 
the anomaly depth and width are mapped into an EFBH 
diameter and depth. If the ratio of an EFBH diameter to 
half-max diameter is close to one for a round anomaly, then 
the anomaly is similar to a flat bottom hole (FBH) or a gap- 15 
ping round delamination. Similarly, if the ratio of equivalent 
uniform gap (EUG) diameter to half-max diameter is close to 
one for a round anomaly, then the anomaly is similar to a 
round void with a uniform gap thickness. Using the IR Con- 
trast method, the EFBH depth is an estimation of the anomaly 20 
depth if the anomaly is similar to a gapping delamination or a 
FBH. Similarly, the EUG depth is an estimation of the 
anomaly depth if the anomaly is similar to a void with a 
uniform gap thickness. The IR Contrast methods allow char- 
acterization of an anomaly in terms of its depth, width, and 25 
gap thickness. The IR Contrast software may be programmed 
for many materials and therefore provides a quick approxi- 
mate prediction of the contrast on a given delamination-like 
anomaly. The software can be used as a training tool for the 
application of flash thermography or other non-destructive 30 
evaluation systems. 

The computer-implemented IR Contrast software com- 
prises two methodologies to analyze the flash thermography 
contrast data for characterization of delaminating-like 
anomalies. The first methodology uses measurements such as 35 
a peak contrast, persistence energy time, withholding time, 
and peak time, which can be plotted against the location of the 
measurement. These measurements show some correlation to 
the anomaly depth, width, and gap thickness. The second 
methodology comprises matching the calibration-based pre- 40 
dieted contrast evolution to the measured contrast evolution. 
The matching of the simulated contrast evolution to the mea- 
sured contrast provides size and depth of the equivalent flat 
bottom hole (EFBH) or equivalent uniform gap (EUG) with 
known attenuation. The EFBH (or EUG) depth and diameter 45 
may be correlated to the width and depth of the measured 
anomaly. Comparison of the EFBH (or EUG) evolutions with 
measured contrast evolutions, so as to provide evaluation of 
delaminations in terms of depth, width, and presence of sur- 
face cracks, confirms this correlation. 50 

Another method known as half-max edge detection 
complements the IR Contrast EFBH (or EUG) width predic- 
tions, and the reverse is also true. The difference in the two 
width measurements is considered to be due to non-unifor- 
mity of the gap thickness, value of gap thickness, and length 55 
effect as well as the limitations of the two estimations. The 
embodiments described herein provide information on the IR 
Contrast/half-max analysis of the flash thermography data. 
Comparison of the adjusted EUG diameter and EFBH diam- 
eter with the half-max width allows assessment of the gap 60 
thickness. Comparison of the peak product time with the 
square of the half-max width also allows assessment of the 
gap thickness. 

Exemplary embodiments will now be described with ref- 
erence to the accompanying figures. Like elements or com- 65 
ponents in the figures are denoted with the same reference 
characters for consistency. 


8 

FIG. 1 is a schematic diagram of an exemplary single-sided 
flash thermography system in accordance with embodiments 
described herein. As shown in FIG. 1 , the equipment for an IR 
flash thermography system 105 comprises a flash lamp 
(source of light/heat) 110 , a flash-hood 115 , flash power 
supply/trigger unit 120 , flash duration controller 125 , an IR 
camera 130 for capturing video images, camera data acqui- 
sition electronics 135 , and a computer 140 . The computer 140 
is used for controlling the flash trigger 120 for the flash lamp 
110 , data acquisition of the camera video data, data display, 
and post-processing of the acquired data. The flash-hood 115 
may be made from sheet metal to form a box-like structure, 
although other structures that form an interior volume, such 
as a cylinder, hemi-sphere, etc. may be used. However the 
flash-hood is structurally configured, one portion of the hood 
may have a large opening so as to enable the hood to be 
positioned on a surface supporting a test object 145 . The side 
opposite to the large opening of the hood has a hole in the 
center to provide a window for the lens of the IR camera 130 , 
which is positioned outside of the flash-hood 115 . The IR 
camera 130 is focused at the top surface of the test object 145 
(part) positioned inside the hood. At least one flash lamp 110 
(two are shown in FIG. 1) is located within the inner wall of 
the flash-hood 115 in proximity to the IR camera 130 . The 
flash lamp(s) 110 directs illumination towards the location of 
the test object 145 without directly shining light upon or into 
the IR camera 130 . The flash-hood 115 functions as a housing 
that contains most of the intense flash. 

If the size of the test object 145 can be accommodated 
inside the hood 115 , then the test object 145 is positioned near 
the hood opening slightly inside the interior volume of the 
hood 115 . Otherwise, the test object 145 is located slightly 
outside of the hood 115 opening. A short duration (e.g., 5 
msec), intense (12 kJ) flash is triggered by the computer 140 . 
The data acquisition is triggered a few seconds before the 
flash and it continues until the prescribed time. The camera 
130 provides a sequence of IR images (or frames) called the 
data-cube (or digital video) of the test object’s surface taken 
at the chosen frame rate (e.g. 60 Hz or 60 frames per sec). The 
intensity (numerical value) of each pixel in the image is a 
function of the surface temperature of the corresponding area 
on the test object 145 at the time of the image frame. The flash 
causes the surface of the test object 145 to warm up slightly, 
and the heat starts to dissipate rapidly. The surface cools 
through thermal radiation, convection, and conduction. Heat 
conduction is considered to be the dominant heat transfer 
mode within the test object 145 until the temperature gradi- 
ents within the test object become small. At later times, the 
heat conduction is of the comparable magnitude of order as 
the combined effect of heat convection and radiation. The IR 
data acquisition and data analysis utilizes the thermal data in 
the short duration immediately after the flash wherein heat 
conduction dominates the thermal dissipation within the test 
object 145 . 

The heat exchange across the boundaries due to the con- 
vection can be assumed to be zero if the Biot number 
(Ng ; =hL/k) is less than 0.1. For example, a half centimeter 
thick graphite/epoxy (k=0.64 W/mK) plate, using h=10 
W/m 2 K, has a Biot number of 0.078. Therefore, heat conduc- 
tion is the dominant mode of heat transfer for this graphite/ 
epoxy plate example. Thinner test objects tend to equalize the 
temperature within the test object quickly and have a rela- 
tively longer cooling time by heat loss to environment. Such 
test objects are not suitable for the IR Contrast simulation 
analysis. A term called threshold thickness (d threshold ) may be 
used to classify a test object as too thin or too thick for heat 
conduction in flash thermography. 
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FIG. 2(a) is a schematic diagram of a plate serving as the 
test object 145. FIG. 2(a) depicts a gapping delamination or 
anomaly 220 with measurement and reference regions of 
interest (Rills) in accordance with embodiments described 
herein, whereas FIG. 2(b) shows a single infrared (IR) video 
image of a test object with a flat bottom hole in accordance 
with at least one embodiment. Referring now to FIGS. 1 and 
2(a), the test object 145 is made of a thermally isotropic 
material with constant thickness 205 and fits inside hood 115. 
The test object 145 is supported at the corners on insulating 
standoffs 150 and the hood is oriented vertically thereabove. 
Assuming that the flash intensity is uniform over the test 
object’s top surface 210, the heat conduction will be in a 
direction normal to the top surface 210 in most of the acreage 
area (i.e., area away from edges of the test object and flash 
boundary). The heat is conducted uniformly from the top 
surface 210 to the bottom surface 215 of the test object. 
Normal heat conduction will be obstructed by an anomaly 
such as a small round gapping delamination 220 at the center 
of the test object, as depicted in FIG. 2(a). The delamination 
220 is characterized by noting it begins at some depth (d) 221 
below the top surface 210. The delamination 220 may also be 
characterized as having some diameter (D) or width 222 and 
a gap 223. The volume bounded by the anomaly 220 on one 
side and the top surface 210 on the other side is called the heat 
trapping volume 225. 

The top surface area 230 surrounding the anomaly 220 
cools faster than the top surface area (or footprint) 235 above 
the anomaly 220. The IR camera 130 captures the surface 
temperature image in terms of the pixel intensity and shows 
the anomaly 220 as a hot spot (e.g. an area warmer than the 
surrounding area) which is about the size and shape of the 
anomaly footprint. The relative pixel intensity of the hot spot 
changes with the time. Deeper anomalies appear at later times 
in the IR video compared to the near surface anomalies. After 
the appearance of an anomaly in the IR video, its relative pixel 
intensity continues to increase with time. The relative pixel 
intensity associated with the anomaly reaches a peak at a 
certain time, and then the relative intensity decays until the 
indication area temperature and the surrounding area tem- 
perature become equal. FIG. 2(b) shows one of the frames in 
a sample IR data-cube from a reinforced carbon-carbon 
(RCC) part with round FBHs machined front the back side to 
simulate gapping delaminations of differing depths and diam- 
eter. IR indications 240 from different flat bottom holes are 
clearly visible in the IR image. A measurement ROI 245 is 
shown as a square in the middle of one of the anomalies. A 
reference ROI 250 is shown as a square outside the anomaly 
area. The part continues to cool down to ambient temperature 
through heat convection and radiation. If any objects contact 
the part or the part has significant area that does not receive 
the flash, then the heat conduction to the unflashed areas or to 
the contacting object also causes cooling. 

In some embodiments, such as for when the test object is a 
laminated thick plate with a large diameter anomaly at a 
shallow depth (i.e. the diameter to depth ratio is very large) 
and with the center point of the footprint of the anomaly 
which is on the top or flash illumination surface, the tempera- 
ture decay after the flash at the center point does not depend 
on the diameter of the anomaly. The temperature decay is, 
however, dependent on the depth of the anomaly, especially 
after the time when a heat wave-front reaches the anomaly. 
The time at which the effect of the heat wave-front, reaching 
the bottom surface, is detectable on the top surface is called 
the normal transit time. For the example of an instantaneous 
pulse of heat radiation being applied to the plate, in which no 
heat transfer from the plate to the environment after receiving 


the flash is assumed, then the surface temperature is governed 
by heat conduction within the plate. The prior publication of 
H. S. Carslaw and J. C. Jaeger (Conduction of Heat in Solids , 
Clarendon Press, Oxford, 1959) provides many analytical 
5 solutions for heat conduction problems. An approximate 
solution for the surface temperature rise at the center point as 
given by J. W. M. Spicer (Active Thermography for Manu- 
facturing and Process Control, Short Course Notes SCO 5, 
Thermosense, SPIE — The International Society for Optical 
Engineering, April 1995) is, 
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where AT is the temperature change from the pre-flash tem- 
perature on the front surface (° C., K, or 0 F.), t is the time (sec) 
measured from end of the flash, x„ is the multiple of heat 
transit time to the anomaly, a is the thermal diffusivity of the 
material or effective diffusivity (cm 2 -sec -1 , m 2 -sec -1 or in 2 - 
sec -1 ), |) is the thermal effusivity of (layered) material above 
35 the anomaly (cal-cm 2 -° C. -1 -sec -1/2 , kcal-m -2 -K. -1 -sec or 
BTU-fT 2 - 0 F. -1 -hr -1/2 ), i s the thermal effusivity of (lay- 
ered) material below the anomaly (cal-cm -2 - 0 C. - -sec - , 
kcal-m -2 -K - 1 -sec - 1/2 or BTU-ft -2 - 0 F. -1 -hr -1/2 ), c is the spe- 
cific heat (cal-gnT 1 - 0 C. -1 , kcal-kg -1 -K -1 or BTU-lbm -1 - 0 
40 F. -1 ), d is the thickness of the (layered) material above the 
anomaly or anomaly depth (cm, m or in), e is the emissivity of 
the flashed surface, At is the flash time (see), H is the average 
heat flux incident on the surface of the test object (cal-cm -2 
sec -1 , kcal-rn -2 sec -1 or BTU-ft -2 hr -1 ), Q is the total heat 
45 incident on the test object per unit surface area (cal-cm -2 , 
kcal-rn -2 or BTU-ft -2 ), k is the thermal conductivity, (cal- 
cm -1 - 0 C. -1 sec -1 , kcal-m -1 -K -1 sec -1 or BTU-ft -1 - 0 F. -1 
hr -1 ), p is the (layered) material density (g-cm -3 kg-m -3 or 
lbm-ff -3 ), and T is the thermal mismatch factor at the anomaly 
50 (T is the Greek letter capital gamma). Equation (1) assumes 
that the surface temperature is affected solely by the diffusion 
of heat within the test object by conduction. The heat con- 
duction at the center point is assumed to be in a direction 
normal to the top surface. The flash from each of the flash 
55 lamps is assumed to approximate the instantaneous heat 
pulse. The anomaly is modeled as a contact between two 
dissimilar layers with a mismatch factor T. For the situation 
with no anomaly and no thermal mismatch, T is equal to zero 
and there is perfect contact between layers of the same mate- 
60 rial. For a transmissive anomaly, T is equal to one and the 
bottom layer is a perfect conductor or heat sink and there is 
perfect contact with the conductor. For the case with a gap- 
ping delamination and no heat transmissivity, the bottom 
layer is a perfect insulator and T is equal to minu s one. It is 
65 assumed that air is present in the delamination gap causing 
heat transmission by conduction across the gap. For thin gap 
delaminations, the delarni- 
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nation gap thickness also factors in the thermal mismatch 
with values ranging typically between minus one and zero. 
Similarly, for bridging within a delamination and uneven 
contact pressure within the delamination, the thermal mis- 
match values would range from minus one to zero, indicating 5 
partial transmission of the heat. For a thick test object, with no 
anomaly, T is equal to zero and equation (1) simplifies to, 


AT(0 = 



(7a) 1° 


Equation (7a) implies that the test object continues to cool to 
its initial temperature indefinitely. The above equations are 15 
not applicable at the instant flash time and earlier times. The 
surface temperature at the “zero time”, which occurs imme- 
diately after the flash, is given by, 


A7o = 



20 

(7b) 


It is difficult to measure the temperature at “zero time” 25 
because it changes very rapidly within the camera measure- 
ment period or the instalment response time. For instance, a 
variation in the delay between flash trigger with respect to the 
camera frame trigger may cause the camera to measure vary- 
ing flash responses. In other instances, some cameras may 30 
saturate momentarily due to the high intensity of the flash. 
Also, the tail end of the flash, if not decayed sufficiently, may 
continue in the first or earlier post-flash frames. Therefore, 
measurements taken immediately after the flash (i.e., the first 
frame) may have lower repeatability and undesired informa- 
tion. Flence, options to block or omit these frames in the IR 
Contrast model and data analysis are provided. For the situ- 
ation in which a large gapping anomaly T is equal to minus 
one, equation (1) simplifies to: 


center point of the anomaly footprint is influenced by the 
lateral (transverse) transit time given by equation (12): 

(12) 

The lateral transit time x h is given by equation (13): 
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where r is the radius of the anomaly (cm, ft) and a l is the 
lateral (or transverse) diffusivity which is the same as the 
thickness diffusivity for isotropic materials . With the IRC on- 
trast model, a single value for the diffusivity that provides an 
acceptable calibration match will be used. This value is likely 
to be between the two diffusivities, the normal and the lateral, 
but can also be dependent on the depth of the anomaly. The IR 
Contrast model addresses the difference in diffusivities 
through calibration constants and limits on the applicability 
of the model. 1 3 is similar to l x and is used conceptually to 
illustrate that the temperature of the center point is affected by 
the diffusion transit times to the lower boundary defined by 
the flaw depth and the edge boundary defined by the radius. 

The effect of the anomaly edge (diameter) on the tempera- 
ture at the center point will occur before the effect of the depth 
of the anomaly ifx, is less thanx. The ratio of these two transit 
times plays an important role in the relative contribution of 
the depth, d, and diameter, D, in the peak contrast time of the 
anomaly. The ratio of the transit times is proportional to the 
square of the diameter to depth ratio, i.e. x,/xa(D/d) 2 . Three 
ratios, gamma (y), omega (to), and phi (i|)), based on the 
dimensional measurements can be defined. Gamma is the 
ratio of the diameter of the anomaly (or the width of a long 
anomaly (j>a3, i|) is defined in Equation (14c)) to its depth 
below the surface as shown in equation (14a), 


D D ap 
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In the above equation, the subscript^ indicates an apparent 
value referring to a value provided by the analysis rather than 
the actual measurement. For an anomaly described by its 
45 depth and diameter, the heat trapping volume is given by 
In some embodiments, the anomaly appears in the IR ^JiD /4. Omega, the ratio of the anomaly depth to test object 

image video sequence at a time related to the normal transit thickness ratio, is shown in equation (14b), 

time given by equation (9), 
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10 j where d* is the test object thickness. Phi, the aspect ratio of 
the anomaly length (major axis) to the anomaly width (minor 
axis), is shown in equation (14c), 


The normal transit time, x, is related to the Fourier Modulus. 

The intensity of the anomaly indication is strongest at another _ / (14c) 

time called the “optimal time” (or the peak contrast time as 6Q w’ 

defined later) and is given by equation (1 1), 


y.v,/. (x, (it) 

Hie anomaly has a lateral heat flow around its center point. 
The heat is conducted over the edges of the anomaly to cooler 65 
surrounding areas and material below the anomaly. The effect 
of the boundary of the anomaly on the temperature of the 


where w is the width of the anomaly. For a linear anomaly, the 
width is defined as the shorter of the two dimensions along the 
orthogonal axes. Hie length of the anomaly, 1, is defined as the 
longer of the two dimensions along the orthogonal axes. 

Next, the terms “normalized temperature contrast” and 
“normalized pixel intensity contrast” need to be defined. The 
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normalized temperature contrast at an anomaly of interest above anomaly definition, the measurement area can be con- 
may be defined as given in equations (15a) and (15b), sidered to be a plate with thickness equal to the anomaly depth 

and the reference area can be considered to be a plate with 
infinite thickness. It was assumed that there was no heat loss 


(T-T°)-(T ref -T^) (15a) 

“ (T-TOi + Wr'f-T^f 

(15b) 

CV = x 100, 

(T - T°) + (Tref - 

where, TJ is the normalized temperature contrast or modula- 
tion temperature contrast and is a fimction of time, t. Time is 
measured from the moment of the flash where time equals 
zero seconds. C % * is the normalized temperature contrast 
expressed as percentage. T is the surface temperature at the 
measurement ROI in ° C., K, or ° F. T re yis the surface tem- 
perature at the reference ROI which is usually outside the 
footprint area of the anomaly. T° is the surface temperature of 
the measurement ROI before the flash. T re ° is the surface 
temperature of the reference ROI before flash. Defining the 
temperature rise or normalized temperature at the measure- 
ment ROI at time t as AT=T”=T-T° and the temperature rise 
or normalized temperature at the reference area outside the 
anomaly at time t as AT,.^T,.^/'=T,. cy -T,. c/ 0 , equation (15a) 
can be expressed as shown in equation (16a): 

, = AT-AT„ f (16a) 

at+at^' 

Thus, the “normalized temperature contrast”, as defined 
herein, is a measure of the ratio of the difference between 
temperature rises (i.e., the relative change in temperature) 
occurring at the anomaly measurement ROI and the reference 
ROI compared to the sum of the two temperature rises occur- 
ring at the anomaly measurement ROI and the reference ROI. 
This definition is similar to the definition of signal modula- 
tion, and the values would range from minus one to one in 
normal, customary circumstances. In some abnormal situa- 
tions, values beyond the customary range are possible. In 
most measurement situations, the reference area is always 
cooler than the anomaly area and, therefore, the contrast 
should normally be positive in value. But in practice, surface 
heat gradients caused by factors such as the proximity of the 
anomaly to the edge of the flash zone or the uneven thickness 
of the cross section may cause the anomaly to appear cooler at 
later times compared to the temperature of reference area. 
The pre-flash surface temperature distribution and emissivity 
variation may also affect the contrast. Ideally, T°=T re ° and 
equation (15a) simplifies as shown in equation (16b): 


T-T ref (16b) 

N ~ T + T rlf - 2T° ' 

The denominator in the contrast definition normalizes the 
numerator. The contrast changes with time and increases 
rapidly after the flash. The contrast reaches a peak value and 
then slowly decays and achieves a substantially constant 
value. The contrast should ultimately reach zero when the 
temperature equalizes in the test object; however, the denomi- 
nator also becomes quite small due to the low temperature rise 
increasing the noise in the contrast evolution at times beyond 
the end (analysis) time (which will be defined later). In the 


5 from the test object and no thermal connection between the 
reference area and the anomaly area. Therefore, the contrast 
would not reach a zero value but would mostly level out or 
reach close to a steady state value after a certain time. 

Most radiometric capable IR cameras are calibrated to 
10 measure the surface temperature. These cameras provide a 
surface temperature map using either a color pallet or a gray 
scale. The contrast formula is directly applicable in this situ- 
ation. However, the IR cameras used in a flash thermography 
set-up may not be calibrated for measurement of the tempera - 
1 5 ture. The IR cameras used in flash thermography may have a 
linear pixel intensity response with the irradiance or the total 
radiation incident on the camera detector pixel. In this situa- 
tion, the total radiation incident on the detector has three main 
components: (1) an emitted radiation from the surface of the 
20 test object; (2) a reflected radiation from the surface of the test 
object; and (3) a radiation from the air between the test object 
and the camera. The main component is assumed to be due to 
the thermal emission as a result of the high emissivity of the 
test object. Therefore, a normalized pixel intensity contrast 
25 may be defined as shown in equation (16c): 


AW-AW ref (16c) 

* “ AW +A W ref ' 

30 

wherein AW re/ =W re/ " =W,. ?/ -W rey ? is the rise in the pixel 
intensity or normalized pixel intensity at the reference ROI at 
time t, AW=W”=W-W° is the rise in the pixel intensity or 
35 normalized pixel intensity at the measurement ROI at time t. 
Thus, the “normalized pixel intensity contrast”, as defined 
herein, is a measure of the ratio of the difference in pixel 
intensity rise (i.e. relative pixel intensity rise) between the 
anomaly measurement ROI and the reference ROI to the sum 
40 of the two pixel intensity rises. 

The pixel intensity registered by the camera due to the 
thermal emissions is approximately proportional to the fourth 
power of the surface temperature of the areas of the test obj ect 
imaged in the pixel . However, it is expected that if the contrast 
45 definition is applied to the total incident radiation (or pixel 
intensity), the contrast would be approximately the same as 
the contrast calculated from the true surface temperatures 
(e.g. C‘=„f) for high emissivity (>0.95) surfaces, assuming 
that the thermal emissivity (within the range of camera sens- 
50 ing wavelength) is high (>0.8) and constant over the surface 
of the test object under inspection. In addition, it is assumed 
that the flash intensity on the surface of the test object, the 
flash duration, and the part/camera/hood geometry are kept 
constant. With these assumptions, the pixel contrast and the 
55 temperature contrast are approximately equal to each other. 
The influence of small variations in the emissivity will be 
modeled approximately using a quantity called the emissivity 
factor ratio (to be defined later). 

The flash time marks the beginning of the contrast evolu- 
60 tion, where normally the contrast value is zero. On acreage 
delaminations, the contrast rises rapidly above the zero level 
with time to reach a peak and then slowly decays to the zero 
level at long times. An ideal shape of the contrast evolution is 
similar to a Gaussian distribution that is skewed to the left. 
65 Confidence in the analysis of the contrast evolution is higher 
if the shape of the contrast evolution is close to the ideal 
shape. Factors such as the part curvature, proximity to the part 
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edge, step in the part, and uneven flash intensity may cause the 
contrast evolution to have a non-ideal shape. Some factors 
cause lateral heat flows that are not related to the anomaly 
under investigation causing a non-ideal contrast shape. In 
such situations, one may instead extract the contrast evolution 5 
with the proper choice of measurement and reference ROIs 
such that the shape of the evolution is as close to the ideal 
shape as possible. In the embodiment described herein, a 
method of artificially leveling or shifting the evolution verti- 
cally is provided, but it should be used only after validation. 10 

In some embodiments, the contrast may not start at the zero 
level. The contrast may also dip below the zero level after 
reaching the peak and then slowly come back to the zero level 
as the temperature equalizes. In this situation, there may be a 15 
time period beyond the peak time when it is below the zero 
level but later it levels out providing a steady value for a 
considerable time. The steady state contrast is given by equa- 
tions (17a) and (17b), 
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(17b) 


The superscript S is used to indicate a period of time in which 30 
the change in the contrast over passage of time is small. This 
period of time is called the “steady state time”, where AT' S and 
AT,/ represent the temperature rises at steady state time for 
the measurement ROI and the reference ROI, respectively, 
and AW S and AW,/ correspondingly represent the pixel 
intensity rises at steady state time for the measurement ROI 
and the reference ROI, respectively. To achieve leveling of the 
tail portion of the contrast evolution, an adjustment to the 
contrast may be applied using the expressions in equations 40 
(18a) and (18b): 
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The Fourier Modulus in heat conduction is defined as the 
variable t' and is called the normalized time. It is related to the 
time as expressed in equation (19): 

t'=th. (19) 

The time ratio of equation (1) can then be expressed as equa- 
tion (20) 



r 



t’ 


( 20 ) 


Substituting equation (20) into equation (1) results in the 
equations of (21a) and (21b): 
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The right hand sideof equation (21b) is only a function of the 
normalized time. Equation (7a) can be expressed as shown in 
equation (22), 
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Substituting equation (21a) and equation (21b) into equation 
(16a) yields equation (23): 
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Selecting a value of 1 4 equal to one would result in bringing 
the contrast at the steady state time to zero. Under ideal 
conditions of flash uniformity and test object geometry, the 
contrast on the flat bottom holes is near zero at steady state 
times. In other embodiments, choosing the reference ROI 
closer to the anomaly reduces the negative incursion of the 
contrast. Uniformity of the flash intensity also helps in reduc- 
ing the negative incursion of the contrast. The flash intensity 
is usually most uniform in the center of the flash-hood. Locat- 
ing the anomaly away from edges of the flash boundary also 
helps in mitigating the negative incursion of the contrast. The 
negative incursions of the contrast may be minimized by 
adjusting the physical test parameters before applying any 
contrast adjustment. 


The normalized contrast is solely a function of the normalized 
time. It should always be positive per equation (23). Equation 
(23 ) implies that the normalized contrast evolution against the 
normalized time is the same for all isotropic materials and 
anomaly depths. Therefore, the evolution is called the invari- 
ant contrast evolution. The normalized contrast is also con- 
sidered to be the thermal response of the anomaly. Equation 
(23) is plotted in FIG. 3 in semi-log and linear fomi. As can be 
seen in FIG. 3, there is an inflection point at t' equal to 0.5. 
Defining an early detection time (t appear ) as one percent con- 
trast. then t' equals to 0.2 and 1 1 =0.2 per equation (9). Here, 
the normalized contrast keeps increasing with time and is 
asymptotic to one. The above contrast evolution is applicable 
to high values of the diameter D to depth d ratio or gamma y. 
Essentially, the contrast evolution is drawn from the center 
points of two plates of the same material. The measurement 
ROI is on a plate with thickness equal to the anomaly depth. 
The reference ROI is on a plate with infinite thickness. Due to 
separation of the plates, there is no cross heat flow between 
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the two plates and the two plates cannot reach the same 
temperature after the flash. While the t hinn er plate reaches a 
constant temperature, the thicker slab continues to cool down 
indefinitely. The contrast evolution for any isotropic material 
with a large FBH-like anomaly can be obtained by using the 5 
above contrast evolution and by transforming the normalized 
time to the real time by using equation (19). The concept of 
using the invariant contrast evolution for prediction of the 
normalized contrast is developed further to design the IR 
Contrast prediction model. to 

As the y ratio becomes smaller, the lateral heat conduction 
contributes to the cooling rate and the contrast evolution has 
a shorter peak. In some embodiments, the contrast rarely 
approaches the value of one, because heat from the trapped 
volume is dissipated to the surrounding area due to cross heat 15 
flow. For delaminations, heat from the trapped volume is also 
conducted to material below the delamination through the 
delamination gap. The peak in the contrast is reached below a 
contrast value of one with a subsequent gradual reduction in 
the contrast to zero due to temperature equalization by the 20 
cross anomaly heat flow. At later normalized times, the heat 
loss mechanisms influence the surface temperature as much 
as the internal conduction. In most embodiments, the contrast 
rarely goes above 0.7 on true delaminations, particularly if the 
reference area chosen is significantly larger (10 times) than 25 
the size of texture and the measurement area is larger (>2 
times) than the texture size. 

In order to accommodate lower values of the y ratio and to 
apply the contrast evolution prediction to practical situations, 
a measured contrast evolution curve (normalized contrast 30 
versus time) taken on a flat bottom hole (FBH) in an isotropic 
material (e.g., an anomaly with a 1 inch diameter and a depth 
of 0.05 inch, resulting in a y reJ ,=20) is used to generate a 
“reference contrast evolution” curve to capture effects of 
most factors in the contrast measurement. The reference con- 35 
trast evolution is considered to be a modified instance of the 
invariant contrast evolution of equation (23). The reference 
contrast evolution curve would be applicable to a range of y 
ratios in the neighborhood of its reference y value. FIG. 4 
shows the normalized reference contrast evolution curve. The 40 
reference contrast evolution curve is similar to the invariant 
contrast evolution curve of FIG. 3 in that it is the source to 
generate the simulated contrast evolution curve, but the pro- 
cess also uses calibration and empirical relationships, which 
will be explained in further detail later. 45 

In the above example, 1 7 equals 0.2 and 1 2 equals 0.35, per 
calculations using equation (9) and equation (1 1 ), respec- 
tively. All material properties as well as the surface emissivity 
are known. Transforming the actual time to the normalized 
time generates the normalized contrast versus the normalized 50 
time curve. In order to generate the simulated (predicted) 
contrast evolution of a flat bottom hole, the material proper- 
ties of flat bottom hole depth and diameter need to be known. 
The normalized time (with known contrast) is transformed to 
the real time. Thus, the contrast versus the real time evolution 55 
is obtained. Making numerous changes to this approach 
allows calibration on known contrast evolutions. Here, the 
term “calibration” refers to the process of defining values of 
certain calibration constants to obtain a better match between 
evolutions of the predicted contrast and the measured con- 60 
trast. 

Using the reference normalized contrast (contrast with nor- 
malized time) and the ratio y, and the three calibration con- 
stants (which will be described and defined later), the pre- 
dicted contrast evolution may be fine-tuned within the 65 
diameter and depth range that is validated during calibration. 

To accommodate the smaller diameter of the anomaly, the 


normalized time expression in equation (1 9) is further modi- 
fied as expressed in equation (24a): 


C'i ( e y ) / a e t • 
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(24a) 


where a e is the effective diffusivity (cm 2 sec x ). The trans- 
formation coefficient Cyis defined as given in equation (24b), 


C f = 


C 3 ^) 
(t 


(24b) 


and used to simplify the expression of the normalized time 
expression as shown in equation (24c), 


f = + C 7 j + C 8 . 


(24c) 


To set-up the calibration based on the diameter of an anomaly, 
a diameter constant C D is introduced as expressed in equa- 
tions (25a) and (25b), 

Y'=C.dY (25 a) 

and 


D^CdD, (25b) 

where C D is the diameter or width calibration constant and the 
nominal value for the diameter constant is equal to one. The 
contrast evolution shape (orprofile) constants, C 3 , C 6 , C 7 , and 
C 8 are derived empirically. Here, “ a e ” is defined as the effec- 
tive diffusivity or an estimate of the effective diffusivity rather 
than the actual diffusivity. The effective diffusivity refers to a 
value of diffusivity used in this model that provides an accept- 
able prediction match after calibration. For isotropic materi- 
als, the effective diffusivity is the same as the material diffu- 
sivity. For anisotropic materials, it is related to material 
diffusivity and to direction of the heat flow. An appropriate 
value of the effective diffusivity may provide adequate cali- 
bration on selected flat bottom holes of a test object. “C d ” is 
the depth (or diffusivity) calibration constant. It is assumed to 
be constant for a given isotropic material within the validated 
range of flaws. The nominal value for the depth constant is 
equal to one. An expression for the simulated (predicted) 
contrast evolution is shown in equations (26a), (26b), (27a), 
and (27b), 


C„J=VC A C 5 A c C^ re }”f 

(26a) 


(26b) 
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where C 5 is the contrast evolution shape (or profile) constant 
and is derived empirically, e" is the emissivity factor ratio 
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used to account for differences in the emissivity between the 
calibration standard and the test object, e' is a ratio of the 
normalized pixel intensity contrast to the normalized tem- 
perature contrast. Higher emissivity or flash power provides 
higher temperature rise and better signal-to-noise ratio. 5 
Lower emissivity or flash power provides lower temperature 
rise and lower value of the simple contrast. Lower emissivity 
also reflects more heat. The reflected heat causes the image 
contrast to decrease due to the increase in the denominator of 
equations (16a), (16b), and (16c). The emissivity factor ratio to 
is introduced: to approximately account for small change in 
the emissivity or power. The value of the emissivity factor is 
equal to one for the object emissivity of one. Usually, the flash 
power setting is not changed and the emissivity of the test 
article is required to be high (>0. 8) and the same as that of the 1 5 
calibration standard. Under these circumstances, the emissiv- 
ity factor ratio is assumed to be unchanged with a value of 
one. The emissivity factor ratio is not needed for prediction of 
the temperature contrast, provided that the part emissivity is 
assumed to be uniform. 20 

Starting with the realtime t (e.g. a frame at 14o or I /120 sec), 
the realtime is transformed to the normalized time t' using 
equation (24b) and then t' is “inverse transformed” to the 
realtime t re/ for a reference material flat bottom hole. Next, the 
contrast is read from the reference curve of FIG. 4 . The 25 
reference curve is stored as a mathematical function F)!,.^). In 
some embodiments, it is stored as a digital curve. The math- 
ematical function when matched is considered to be a curve- 
fit and allows for easy calculation of derivatives (e.g., F'(t,. e/ -), 
F"(X re y), etc.). The transformation/calculation process steps 30 
are given below. The reference contrast evolution may be 
obtained from a material different from the non-destructive 
evaluation (NDE) reference standard. The transformation is 
indicated by— > and the calculation is indicated by => . The 
realtime transformation is expressed in equation 28, 35 

(28) 

A flaw size parameter (FSP), A c , defined as a function of y, is 
expressed in equation (29): 


20 

constant calibration curves established from flash thermog- 
raphy data on a calibration test sample having flat bottom 
holes. C 5 is a contrast evolution shape constant based on 
values of the contrast evolution shape constants C 4 and C 9 . C 5 
has a normalizing effect of the FSP and is expressed as 


1 ( 30 ) 


In equation (30), the following values were substituted, 
d=0.05 inch, D=1 .0 inch; and y'=D/d 20. Values of the evolu- 
tion shape constants are derived empirically. A flowchart of 
the process of generating the calibration curves is provided 
and described in further detail with reference to FIG. 30. A 
flow chart exemplifying the steps of the method for the above 
contrast simulation is shown in FIG. 10. The method begins at 
step or block 1 005 , where material properties, depth constant, 
anomaly properties, frame rate, and the calibration file are 
used as inputs. At block 1010 , the flaw size parameter A^ is 
calculated by using equation (29). Once A c is calculated, its 
value is used to derive values of amplitude constant and 
diameter constant from the respective calibration curves. At 
step or block 1015 , an array of times (t) at the chosen frame 
rate and desired number of frames is made and then trans- 
formed to f from equations (24a) and (24c). At step or block 
1020 , the hardcoded normalized reference contrast evolution 
is used to read the reference contrast for each f as given by 
equation (28). At step or block 1025 , the contrast for each t is 
computed using equations 26(a) and 26(b). At step or block 
1035 , the simulated contrast evolution is displayed. 

In one embodiment, a direct depth estimation can be made. 
By manipulating equation (24a), an expression for the depth 
of an anomaly can be obtained as expressed in equation (31), 


Jx'G) (29) 

Ac = 1 -e U 4 J, 0<A C < 1. 

C 4 and C 9 are the contrast evolution shape constants. The FSP 
(also called “parameter 1”) is used for amplitude and diam- 
eter calibration. FIG. 5 shows a plot of the FSP as a function 
of the depth and diameter of the anomaly. FIG. 6 shows a plot 
of the FSP as a function of ganuna. For high values of gamma, 50 
the value of the FSP is equal to one. FIG. 7 shows the trans- 
formation coefficient C^-as a function of the flaw size param- 
eter A c . For high values of the FSP, the value of the transfor- 
mation coefficient C f is equal to one. The value of the 
transformation coefficient increases with a decrease in the 55 
value of the FSP. The inverse of the transformation coefficient 
is plotted against gamma in FIG. 8. This plot indicates that, 
for high values of gamma, the value of the transformation 
coefficient is equal to one. 

During the calibration process, the amplitude constant (C A ) 60 
and the diameter constant ( C D ) are determined by trial and 
error for a given flat bottom hole. Many flat bottom holes are 
used to establish the calibration/correlation relationship 
between the FSP and the amplitude constant and between the 
FSP and the diameter constant. Usually a curve fit in the 65 
calibration data is used to establish these functions. FIG. 9(a) 
is a representative plot of the diameter constant and amplitude 


, t -U (3D 

J_ = (f'-C 8 )>l- e c >) 1 
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C 3 \er) 

By estimating y’, equation (31) can be used at many t and f 
values to estimate the depth. Starting with the measured con- 
trast data, which would typically be a two-column array (col- 
umn 1 provides the time or the frame number, column 2 
provides the contrast) with a time spacing of 0.01 7 sec (Vfco of 
sec based on 60 Hz frame rate), the peak contrast of the 
measured data is found. Using equation (26a), a multiplier, 
M, can be defined as expressed in equation (32), 

q t-p*°k g-pcak ( 32 ) 

M s s' £C a C 3 A c = *° , = — . 

h -fS-ptak 0.56 

'w.rrf 

For equation (32), the subscript “a” indicates the actual or 
measured contrast. The denominator with a value of 0.56 
comes from the peak amplitude of the invariant evolution 
derived from FIG. 4 . By manipulating equation (32), A c can 
be expressed as shown in equation (33), 
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A c = . 

£'?C„C 5 

5 

By manipulating equation (29), y' can be expressed as shown 
in equation (34), 

_L (34) 

r' = (-(ln(l-4 c ))C 4 ) c 9. 10 

By manipulating equation (25), the diameter can be expressed 
as shown equation (35), 

15 

n - — (35) 

“ <V 
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Equation (35) can be used to estimate the equivalent diameter 
of the flat bottom hole, if all quantities on the right hand side 
are known. Equation (34) cannot be used to estimate y’ until 
we know the value of A c . Per equation (33), the value of A c is 
not known until the value of C A is known. However, an itera- 2 s 
tive loop can be set up between equations (31) through (35). 
The starting values for C D and C A are the nominal values 
equal to one. M is calculated using equation (32) as a constant 
for the iterative loop. Equation (33), followed by equation 
(34), is used to provide the value of y'. At this point, the 30 
calibration constant for C D and C A are calculated, as they are 
functions of A c . Next, equation (31) is used to estimate the 
depth, and equation (35) is used to estimate the diameter. 
Next, A c can again be estimated using equation (33) with the 
latest value for the amplitude constant C A , and the process is 35 
repeated until the change in the depth estimation in successive 
iterations is less than a predetermined set value. The flowchart 
of this reiterative process is illustrated in detail in FIG. 11 , and 
the process is called the direct estimation of the equivalent 
FBH depth and diameter. The method begins at step or block 40 
1130 , where the contrast evolution is computed, followed by 
filtering the contrast evolution at block 1135 . At step or block 
1140 , the peak amplitude of the contrast evolution is mea- 
sured and the multiplier M is calculated using equation (32). 

At step or block 1145 , setting the value of C A equal to one for 45 
the first iteration, the flaw size parameter A c is calculated 
using equation (33). At step or block 1150 , the adjusted diam- 
eter to depth ratio is calculated using equation (34). At step 
1155 , setting the value of C 4 equal to one for the first iteration 
and using the value of A c from equation (33), the anomaly 50 
depth is calculated. At step 1160 , the anomaly diameter is 
estimated using equation (35). At decision block 1165 , if the 
change in the depth estimate from the previous estimate is less 
than a predetermined set value, then the direct depth and 
diameter values are finalized. If not, the method loops back to 55 
step 1145 andre-perfonns the calculations through step 1160 
until the change in the depth estimate is less than the prede- 
termined value specified at step 1165 . 

The depth estimation may be done at numerous times dur- 
ing the measured contrast evolution to improve confidence in 60 
the value derived. For instance, it can be done at the peak 
contrast time as shown for Calculation number 1 in the table 
of FIG. 12 . The other points of interest include the left 50 
percent peak contrast point (Calculation number 2 values of 
FIG. 12 ), the right 50 percent peak contrast point (Calculation 65 
number 3 values of FIG. 12 ), and the right 10 percent peak 
contrast point (Calculation number 4 values of FIG. 12 ). For 
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time beyond the right 50 percent contrast point (t 5aR ), influ- 
ence of the lateral heat diffusion processes may be noticeable 
as a significant change in the contrast evolution. 

The IR Contrast model uses several constants in equation 
(24a) and equation (25a). Equation (24a) is the time transfor- 
mation equation. Therefore, the constants appearing in equa- 
tion (24a) are called the “transformation constants”. These 
include C 1; C 3 , C 6 , C 7 , C 8 , C d , and C n . The constants C 15 C 3 , 
C 6 , C 7 , and C 8 have fixed values . Constants C 7 and C 8 are used 
to set the time offset in the transformation. Constants C 1: C 3 , 
and C 6 are used to manipulate effect of y. Hie effect of change 
in y is not fully accounted for by these internal constants, 
especially for the lower y values. The effect of change in y on 
the shape of the contrast evolution curve (e.g. the peak ampli- 
tude, the peak time, etc.) may need adjustment or calibration 
for change in the diffusivity from material to material due to 
change in anisotropy and layer thickness. Therefore, an exter- 
nal constant C D as a multiplier to y is introduced. This external 
constant is calibrated to capture the behavior of measured 
contrast on a given material, and therefore it is considered as 
one of the calibration constants. Since C D multiplies to y, it is 
also called the “diameter constant”. The effect of change in 
C D is more complex. It changes both the amplitude and the 
peak time (or time scale) of the contrast evolution and 
changes the estimation of the EFBH diameter. 

The external constant C d is used as a multiplier to the 
diffusivity. Change in C d causes a change in the time scale or 
the shape of the profile but the peak amplitude is not affected. 
The constant is used to choose (or calibrate to) the effective 
diffusivity of the material and can be determined by matching 
the simulated contrast evolutions to the measured contrast 
evolutions of the FBHs with known diameter and depth. It 
multiplies to the depth (1/d 2 ). Therefore, it is also called the 
depth constant. Because C d also multiplies to the diffusivity, 
it can also be called the diffusivity calibration constant. The 
effect of change in the diffusivity in analyzing a given contrast 
data results in a change in the time scale and the depth esti- 
mation of the anomaly. The calibration process is intended to 
minimize errors due to small differences between the actual 
diffusivity and the effective diffusivity. Equations (26a) and 
(26b) involve additional constants. They include C A , A c , and 
C 5 . These quantities have the effect of changing the amplitude 
of the contrast (amplification of vertical scale) and have no 
effect on the time. Therefore, C A and C 5 are called the “ampli- 
tude constants”. Constant C 5 has a fixed value. 

The flaw size parameter is not a constant. A c is a function 
of y . Internal constants C 4 and C 9 are used in the computation 
of A c . The effect of change in y on the contrast amplitude is 
not fully accounted for by A c . Also, the effect of y on the 
amplitude may need adjustment or calibration for change in 
the diffusivity from material to material. Therefore, an exter- 
nal constant C A , used as a multiplier to A c , is introduced. This 
constant is calibrated to capture the behavior of the measured 
contrast amplitude, and therefore it is considered as one of the 
calibration constants. The effect of C A is to change the ampli- 
tude without affecting the time scale. Therefore, C A is called 
the amplitude calibration constant. 

Constants C 4 through C 9 are internal constants with fixed 
values. They are used in the time transformation as well as the 
computation of the contrast amplitude. Internal constants are 
not part of the calibration constants. The IR Contrast uses 
three external constants for calibration. They are C A , C D , and 
C d . The calibration constants capture effects of variation in y 
or the flaw size parameter, A,., on the contrast evolution. 
Therefore, calibration profiles of C A and C D , as functions of 
A c , are established. The nominal value for each of the con- 
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stants is one, and usually for high values of the A c (>0.5), the 
corresponding values of C 4 and C D are expected to be close to 
the nominal value of one. 

A typical calibration process for an isotropic material 
involves use of various acreage FBHs in the subject material 5 
providing the A c values from 0.01 to at least 0.5. An optimized 
contrast evolution is obtained on each of the flat bottom holes. 
Each contrast evolution is analyzed separately starting with 
the highest A c (>0.5) values. The data contrast evolution is 
displayed in the IR Contrast application, and a simulation 10 
contrast evolution is generated by choosing inputs for mate- 
rial, emissivity factor ratio, FBH depth, and FBF1 diameter 
(with depth constant=l, diameter constant=l, and amplitude 
constant. The process of matching simulation reveals any 15 
changes needed in the depth constant. The exercise is 
repeated on contrast data from many high flaw size parameter 
value holes, and a value of depth constant providing the best 
match on the contrast evolution curve is determined. Once the 
depth constant has been determined, contrast evolution 20 
curves with lower flaw size parameter values are analyzed 
with the calibrated (fixed) depth constant. The analysis 
involves matching the simulation with the measured contrast 
by using the actual FBF1 depth, the actual FBF1 diameter, and 
the depth constant as determined above. In addition, the diam- 25 
eter constant and amplitude constant are adjusted in an itera- 
tive manner to provide the best match of the simulation. The 
best match values of the diameter and amplitude calibration 
constants are plotted against the respective flaw size param- 
eter value. Curves fitted through the plotted data are called the 11 1 
calibration curves. The amplitude and the diameter constant 
curves diverge with reduction in the flaw size parameter 
value. 

The emissivity factor ratio e" is normally set to one during , 
calibration instead of choosing the actual value for emissivity. 
Normally, it is impractical to assess emissivity at the mea- 
surement region. It is assumed that the emissivity is high and 
uniform over the surface of the test object. Thus, the emissiv- 
ity factor ratio can be set to one during calibration. Similarly, 40 
the amplitude ratio £ (explained in further detail later herein) 
is set to one during calibration on the FBHs. The emissivity 
factor ratio and the amplitude ratio have identical effect on the 
contrast prediction. The iR Contrast simulation assumes that 
the flaws are located at a shallow depth compared to the part 45 
thickness. It also assumes that the part thickness is considered 
to be quite thick, i.e., much greater than the threshold thick- 
ness. 

Even the approach of averaging four points (to find a better 
match of the simulation contrast with the measured contrast) 50 
is not likely to provide the optimum solution due to differ- 
ences between the shape of the simulated contrast evolution 
and the measured contrast evolution. The direct depth esti- 
mation may provide a solution close to the optimum solution. 
Thus, the direct depth estimation is considered to be a starting 55 
point for a more accurate depth estimation using the follow- 
ing method. 

The second method of estimating depth is based on scan- 
ning possible values of depth and diameter. An array of “1” 
depth values and an array of “J” diameter values is defined. 60 
Both I and J are natural numbers. The direct estimation of 
depth and diameter are captured in these arrays, possibly as 
the center values. This results in I times J diameter and depth 
pairs. The simulation contrast for one depth-diameter pair at 
a time is calculated. The simulation match or the fit index 65 
“delta square” is calculated from the measured contrast using 
the expression of equation (36), 


«/■ 
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A second fit index “delta square average” is calculated using 
the expression of equation (37), 


AC* = - 
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( 37 ) 


A third fit index “percent delta square” is calculated using the 
expression of equation (38), 


A c* * = 100 
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( 38 ) 


where ns and nf are the starting and ending frame numbers, 
respectively. Typically, the range may span from a frame at 
left 50 percent peak contrast to a frame at right 10 percent 
peak contrast. The delta square, equation (36), is calculated 
for all depth-diameter pairs and the minimum of all delta 
squares is sorted out. The diameter-depth input providing the 
minimum delta square is the best possible match within the 
selected range of frames. Equation (37) provides the average 
difference (difference square per frame) as an indicator of the 
fit between the simulation and the measured data. The method 
of array search is called the “scan fit.” When the array is 
constructed around the direct depth estimation solution, then 
the method is termed as the “search fit.” Equation (38) pro- 
vides the percent difference between the simulation and the 
measured data. 

The raw contrast data contains camera pixel noise. The 
signal can be filtered by a signal processing routine to 
improve the IR Contrast calculations. The filtered contrast 
and the measured raw contrast are shown in FIG. 13. Many 
signal smoothing options are available in the IR Contrast 
method. They include a Fast Fourier Transform (FFT) filter, 
an infinite impulse response filter, a finite impulse response 
filter, an exponential running average filter, an alpha beta 
filter, and a Kalman filter. Detecting frames that provide the 
peak contrast, such as left 50 percent of peak frame, etc., are 
affected by the noise incursion in the selected frame. Signal 
smoothing using signal processing filters (e.g., low pass FFT) 
helps to locate the frames more precisely. Also, the delta 
square calculation is affected somewhat by the noise and thus 
preferable to perform the delta square calculation on the 
smoothed (or smooth) data. 

A smoothing index is defined as a sum of square of differ- 
ence between the raw (unfiltered) contrast evolution and the 
filtered contrast evolution using the equation (39), which is 
similar to equation (36). The smoothing index “delta square” 
is calculated using the expression of equation (39), 



25 


US 8,577,120 B1 


nfl 

ACU = Y (C wM -C W J. 


(39) 


5 

A smoothing index “delta square average” is calculated using 
the expression of equation (40), 


a = = 
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A smoothing index “percent delta square” is calculated using 
the expression of equation (41), 
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where nsl and nfl are the starting and ending frame numbers 
used in this calculation. The subscript “af ’ is used to indicate 
the filtered or smooth contrast evolution. Lower values of the 
smoothing index imply a better match but not necessarily 
better smoothing, which is controlled by the degree of 
smoothing input. Starting from a highly low-pass-filtered 
(low-pass frequency f) signal, for any decrease in the filtering, 
(i.e. increase f), the delta square values would correspond- 
ingly be decreased. The derivative of slope, ld 2 AU M , 2 " / /df 2 l, 35 
decreases to zero with the low-pass frequency. The slope IdA 
C w 2/ /dfl reduces with low-pass frequency, and at one point no 
significant reduction in the delta square is observed at the cost 
of increasing the low-pass frequency. Therefore, a compro- 
mise is chosen between the degree of smoothing and the 40 
smoothing index as the optimized smoothing solution. If a 
consistent degree of filtering (filter type and setting) is used, 
then equation (41) can be viewed as a measure of noise-to- 
signal ratio. The measured noise-to-signal ratio is higher than 
this estimation, as the filtered signal is not completely smooth 45 
similar to the predicted evolution. FIG. 13 also shows the 
predicted contrast curve that provides the best match to the 
measured contrast evolution. The EFBH depth and diameter 
input to the simulation providing the best match to measured 
contrast data is obtained in an automated way. The IR Con- 50 
trast method calculates the approximate depth and diameter 
of the EFBH directly based on the peak contrast, peak con- 
trast time and other times including t 50i , t 5aR and t laR . The 
simulation matching provides an indirect method of assessing 
the EFBH depth and diameter. 55 

In accordance with embodiments described herein, a flow- 
chart of an exemplary process of estimating depth and diam- 
eter values via a simulation match is provided with reference 
to FIG. 14. At step or block 1430, the estimation process 
begins by computing the contrast evolution. At step 1435, the 60 
contrast evolution is filtered. Next, either at step 1440A or 
step 1440B, the process comprises the step of assuming the 
initial diameter and depth, or the step of estimating a direct 
depth value, respectively. At step 1445, a range of widths and 
depths for searching about the current diameter and depth 65 
values is established. At step 1450, the simulated contrast 
evolution for all diameter-depth pairs in the range from step 
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1445 are calculated. Next, at step 1455, the delta square 
values between all simulated evolutions and the filtered con- 
trast evolutions are calculated. At step 1460, the simulated 
evolution with least delta square is located. The lowest delta 
square and the corresponding diameter and depth are noted. 
At step 1465, the depth, diameter, and delta square are com- 
pared with the values from the previous iteration. If the 
change is less than a predetermined, specified value, then the 
process continues to step 1470. If not, then the process loops 
back to step 1450 and the step of calculating the simulated 
contrast evolution is repeated along with the steps 1455, 
1460, and 1465 to arrive at compared values again. Once the 
process arrives at step 1470, the final depth and diameter 
values are determined using a simulation match. 

In another embodiment, a separate contrast imaging 
method application is written to work with a commercially 
available IR data analysis software application sold under the 
brand name of EchoTherm®. For this embodiment, the con- 
trast imaging application allows for the choice of the refer- 
ence area, and the method then converts the raw data-cube 
into the contrast data-cube using the reference area chosen. 
The temperature versus time profile of each pixel is replaced 
by the normalized contrast versus time evolution. The con- 
trast data-cube images look similar to the raw data-cube 
images. FIG. 15 shows a frame from a raw data-cube with 
temperature versus time curves at selected regions of interest 
(ROIs), including the reference ROI. Here, all pixels are used 
as the measurement pixels. Thus, the pixels within the refer- 
ence ROI would have a contrast value of zero or close to zero. 

FIG. 16 shows a frame image from a contrast data-cube 
with the contrast versus time curves at selected ROIs. The 
contrast data-cube allows for a quick view of the contrast 
evolution within the vicinity of the chosen reference ROI and 
forms the basis for feature images based on signal processing 
(e.g., filtering) and on the thermal measurements (e.g., peak 
time) described further in the following paragraph. These 
feature values are calculated for each pixel and displayed as a 
feature frame image. In accordance with another embodi- 
ment, the flowchart for the overall. IR Contrast process is 
provided in FIG. 17. The process begins at step 1705, wherein 
the flash thermography video data is recorded as a sequence 
of digital frames. At step 1720, an anomaly on one of the 
digital frames is selected. At step 1725, measurement and 
reference region of interest are determined. At step 1730, the 
contrast evolution data of the anomaly is extracted using 
either a single point extraction for round anomalies or a 
multipoint extraction for linear anomalies. At step 1735, the 
contrast evolution is computed for each measurement point 
followed by the act of filtering the contrast evolution at step 
1740. At step 1745, the peak contrast, the peak contrast time, 
the left fifty percent peak time, and the right fifty percent peak 
time are calculated. From steps 1705, 1725, or 1745, the 
process can continue to step 1760, wherein an appropriate 
frame between the left fifty percent peak time and the peak 
contrast time is located, followed by the actions performed at 
step 1765, wherein the half-max boundary and the half-max 
width of the boundary of the anomaly are calculated. At step 
1745, the process can also proceed to step 1750, wherein the 
peak product time is measured along with other features of the 
contrast evolution, followed by the actions at step 1755, 
wherein the anomaly gap thickness based on correlation is 
assessed by calculating the ratio of the peak product time to 
the square of the half-max width. Alternatively, from step 
1745, the process can proceed to step 1770, wherein an 
equivalent flat bottom hole or equivalent uniform gap diam- 
eter value is estimated, such as by assuming an attenuation of 
zero for the equivalent flat bottom hole or a fifty percent 
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attenuation for the equivalent uniform gap. Then, at step 
1775 , the anomaly gap thickness based on correlation is 
assessed by calculating a ratio of the equivalent flat bottom 
hole diameter or of the equivalent uniform gap diameter to the 
half-max width. Additionally, from step 1705 , the process 5 
may continue to step 1780 , wherein a reference region of 
interest is selected. Continuing with this thread of the process 
at step 1785 , the sequence of digital frames is converted to a 
normalized contrast data-cube. At step 1790 , the features of 
the contrast evolutions for each pixel are measured, which in to 
turn are subsequently displayed as a feature frame image at 
step 1795 . 

The contrast evolution on a FBH has a shape close to the ideal 
shape of the contrast evolution. The ideal shape is similar to a 
skewed statistical distribution. Therefore, measurements 15 
based on the ideal shape of the evolution are possible. The 
shape of the contrast evolution is affected by two pixel inten- 
sity decay curves, e.g. the pixel intensity decay at the mea- 
surement ROI and the pixel intensity decay at the reference 
ROI. FIG. 18 shows the measurement and reference ROI 20 
decay curves plotted in log time. In FIG. 18 , the reference 
pixel intensity versus the log of decay time is more or less 
linear until about frame 248 . After this frame, the slope of the 
curve turns and becomes horizontal. The bend in the slope of 
the curve signifies the point in time of the effect of heat 25 
reaching the lower surface and affecting the heat diffusion 
from the reference ROI. For isotropic materials, the pixel 
intensity log decay is linear until the equilibrium temperature 
(X taiI described in further detail later herein) is reached. In the 
normalized contrast definition, it is assumed that the refer- 30 
ence decay curve with the log of time is linear. The reference 
decay curve is dictated by the heat diffusion or diffusivity in 
the normal direction. Apart from using the simulation con- 
trast, thermal feature measurements on both the measured and 
the IR Contrast simulated evolution can be performed. An 35 
explanation of each measurement parameter now follows: 

Peak Contrast, Cj; peak \ 

The peak contrast is defined as the “amplitude” or value of 
the highest point in the contrast evolution curve. It may also 
be referred to as the peak amplitude or contrast amplitude. 40 
FIG. 19 shows a graph with the peak contrast and peak con- 
trast time identified. The peak contrast is related primarily to 
the dominant feature of the anomaly. When surface texture 
(unlike delamination) is present, it provides its own peaks 
which may be superimposed on the peak caused by the sub- 45 
surface anomaly. The surface texture indications occur at 
earlier times and the subsurface delamination indications 
occur at later times depending upon the depth. The peak in 
contrast is reached at later times if surface texture indications 
are not present. However, the texture indications may provide 50 
significant contrast at earlier times and may advance the 
resultant peak contrast time. Thus, in some embodiments the 
peak contrast is used as a measure of the resultant effect of the 
surface and subsurface anomaly. The peak contrast has a 
strong relation to the anomaly width and gap thickness in the 55 
anomaly. 

Peak Contrast Time (sec), t peak : 

The peak contrast time is defined as the time (flash time=0 
sec) at which the peak in the contrast occurs. The peak con- 
trast time is also called the peak rise time or peak time. The 60 
contrast increases until the peak time and decreases after the 
peak time. Peak time divides the time into two time windows, 
namely the contrast rise time window and the contrast decay 
time window. The time window from the end of flash time to 
the peak time is called the contrast rise time window. The time 65 
window starting from the peak time and thereafter is called 
the contrast decay time window. The peak contrast time is 
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influenced by the superimposition of effects of the surface/ 
near surface texture and sub-surface anomaly. The two mate- 
rial conditions may provide either a single merged peak or 
two distinct peaks in the contrast evolution. It is assumed that 
the peak contrast and peak contrast time represent the domi- 
nant feature of the two effects if the subsurface peak is 
expected to occur at a later time compared to the texture peak. 
Thus, the peak contrast time has a strong relation to the 
anomaly depth. 

Fifty Percent Left of Peak Time (sec), t 50i : 

This parameter is the time (flash time=0 sec) at which the 
50 percent of the peak contrast occurs before the peak contrast 
time. The contrast evolution before t 50£ may be useful if the 
surface/texture anomalies are needed to be analyzed. The IR 
indications of the surface anomalies are sharpest after the 
flash and diffuse away with time. Depending upon the loca- 
tion and size of the measurement ROT with respect to the 
texture indication and the severity of the texture effect, the 
contrast evolution shows variation in early times (<t 50 z) after 
flash. This time has a strongly relation to the anomaly depth. 

Fifty Percent Right of Peak Time (sec), t 50 ^: 

This parameter is a time (flash time=0 sec) at which the 50 
percent of peak contrast occurs after the peak contrast time. 
The contrast evolution beyond this time is significantly influ- 
enced by slow diffusion and cooling processes. The lateral 
heat diffusion, if present, can affect the contrast noticeably 
between the t 5QR and t end . Therefore, a portion of contrast 
evolution between t 50i and t 5aR is recommended to be used 
for comparison with the simulated contrast evolution in 
evaluation of the subsurface delaminations. 

Ten Percent Right of Peak Time (sec), t 10J ;: 

This parameter is a time (flash time=0 sec) at which the 10 
percent of the peak contrast occurs after the peak contrast 
time. 

End Time (sec), X end : 

This parameter is the time corresponding to frame ns in the 
calculation of the delta square. The relationship between the 
end time and fifty percent right of peak time and the peak time 
is given in equation (42a), 

t en<P- t 50R'~ tpeak- (42a) 

Tail Time (sec), t MlZ : 

Another time used in the leveling of the contrast signal or 
in optimizing the location of the reference ROT is the tail 
time. The IR Contrast method assumes that the contrast value 
is nearly zero at the time of t taiI . The tail time is expressed in 
equation (42b), 

^ [ail 3^50 R~^tpeak- (42b) 

Fifty Percent Persistence Time (sec), X 50L _ 50R : 

The time between fifty percent left of peak time and fifty 
percent right of peak time is given by equation (43), 

^50L-50R ^50R~^50L" (43) 

The fifty percent persistence time is strongly related to the 
anomaly width. The times defining the persistence can also be 
based on other time parameters, e.g., t 50i _ laR . 

Withholding Heat Time (sec): 

This parameter is a measure of area under the contrast 
evolution curve and is strongly related to anomaly width. An 
expression is shown in equation (44a), 

t =T' t-peaka 

‘Wj, 50-50 < “w,a l 50L-50R- 


(44a) 
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Persistence Energy Time (sec): 

This parameter is a measure of time given by equation 
(44b) 

5 

"50 r (44b) 

n 50L 

fE, 50-50 - ■ 

V 

10 

where v is the frame rate in frames/sec (Hz). n 5at and n so/f are 
the frame numbers corresponding to the times of t 50i and t saR , 
respectively. The Persistence Energy Time has a strong rela- 
tion to the anomaly width. The frame numbers could be 
different as chosen by the user, e.g^njo^tonjo^njo^ton^^; 15 
or nl to n2. Depending upon the choice of the starting frame 
number and ending frame number in relation to the peak time, 
the integrated times of equation (44b) may be denoted as the 
rise energy time, the decay energy time, etc. In many situa- ., (| 
tions of inspection, especially for thin parts, use of the rise 
energy time may be appropriate, as the rise time is much 
shorter than the decay time and relates to the anomaly char- 
acteristics better due to the higher rate of the temperature 
change. Also, if the anomaly depth is expected to be constant, 2 s 
integrating between fixed frames (e.g., nl and n2) is useful, 
especially within the rise time window. Both the persistence 
energy time and the withholding heat time are related to the 
heat trapping volume, or to the depth and diameter of the 
anomaly. 30 

Two other parameters based on time may be defined. They 
are the “peak product time” and the “peak time to contrast 
ratio” (or peak time per unit peak contrast). The peak product 
time (sec), t pmduct , is given by equation (45a), 

35 

H5 a) 

The peak time to contrast ratio (sec), t r c , is given by equation 
(45b), 

40 

[peak (45b) 

t,C “ j4~ptak ‘ 

L wp 

From the contrast maps shown in FIG. 20, FIG. 21(a), and 45 
FIG. 21(A), it is evident that the peak product time is related to 
the diameter of the EFBH (or EUG) and the peak time to 
contrast ratio is related to the depth of the EFBH (or EUG). 
The contrast map in FIG. 20 is generated by using an image 
post-processing simulation routine of a commercially avail- 50 
able software known as ThennoCalc, sold by Innovation Ltd., 
whereas the contrast maps in FIG. 21(a) and FIG. 21(A) were 
plotted using the calibrated IR contrast of the embodiments 
described herein. FIG. 21(a) is from a flat bottom hole cali- 
bration, and FIG. 21(A) is from the same size (in depth and 55 
diameter) but thins voids with 50% attenuation. Thus, these 
measurements are useful if the simulation analysis is not 
performed. A smoothed contrast evolution and even simula- 
tion fit can provide a more stable detection of the peak point 
in the contrast evolution compared to the peak detection on 60 
the raw contrast evolution. Therefore, for more repeatable and 
reliable measurements, a smoothed evolution or the simula- 
tion fit evolution is recommended. In some circumstances, the 
raw evolution may not have a peak corresponding to the 
subsurface anomaly due to the influence of the texture. In 65 
such a case, use of the simulation fit evolution that is fitted to 
the subsurface portion of the raw contrast evolution may be 
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better due to higher confidence in locating the peak and com- 
puting the quantities in equations (45a) and (45b). 

Adjusted Flaw Size Parameter (adjusted FSP), A cA : 

This parameter is also known as flaw size parameter 2 and 
is given by equation (46), 

(46) 

The adjusted FSP relates to the peak contrast per equations 
(26a) and (26b) and is illustrated in FIG. 22. The data for FIG. 
22 was taken from an RCC sample. The adjusted flaw size 
parameter provides correlation with the peak contrast and 
therefore is useful in statistical analysis, such as probability of 
detection analysis. 

The persistence energy time, withholding heat time, and 
peak contrast are all measures of signal strength. The FSP and 
the adjusted FSP are measures of the anomaly size. Therefore, 
these measurements are candidates for a conventional prob- 
ability of detection (POD) analysis where signal response to 
flaw size data is analyzed. The product of peak time and 
persistence time provides a correlation to the product of depth 
and diameter. The withholding heat time and persistence 
energy time capture the strength of the IR indication to draw 
attention for further analysis including the half-max width 
and the equivalent flat bottom hole analysis. Shallow anoma- 
lies have shorter persistence but higher amplitude. Deeper 
anomalies have lower amplitude. Wider anomalies have 
longer persistence. The persistence energy time captures the 
portion of evolution containing the peak of the contrast. 

In a contrast evolution, the contrast is zero until time 
denoted by t appear . The rise time (which is equal to t t - 
'■appear) is much shorter than the decay time (i l0R -l peak ). An 
inflection point exists between the t peak and t appear . Also, an 
inflection point is located between t peak and t 10J! . Different 
temperature gradients and heat flows exist in different time 
zones of the contrast evolutions of the FBH. Consider a FBH 
with a high y value. For a post-flash time less than t 50i , the 
normal heat flow dominates and is therefore related more to 
the depth of the anomaly. The nomial heat flow is highest 
immediately after the flash and reduces rapidly until the time 
t 50i . The lateral heat flow is negligible immediately after the 
flash till the early detection time ( t appear )• Therefore, the early 
detection time strongly relates to the depth but not to the 
diameter. Unfortunately, it is a bit cumbersome to locate the 
early detection time due to the low contrast amplitude. 

As the normal heat flow reduces, the lateral heat flow 
increases in the trapped volume with time. Thus, direction of 
the heat flow changes from vertical to a near lateral in the heat 
trapping volume. The lateral heat flow causes “spilling” of the 
heat from the boundary of the heat trapping volume into the 
surrounding volume. The normal conduction is dominant in 
the surrounding area for the period of the significant contrast. 
A small value for oo is assumed so that the heat transit time to 
the back surface is much longer than the heat transit time to 
the anomaly or the lateral transit times. Between the time t 50i 
and the \ peak> the normal heat flow reduces significantly and 
the lateral heat flow builds up in the trapped volume. The time 
t 50 L > s close to the inflection point in the contrast evolution, a 
time at which the normal temperature gradient in the heat 
trapping volume becomes low and the contrast continues to 
rise due to outpacing of the cooling in the surrounding area 
compared to the lateral heat flow (spilling). The inflection 
point near the t 50i signifies transition from dominant normal 
heat flow in the trapped volume to the dominant lateral heat 
flow and the resulting heat flow across the boundary of the 
heat trapping volume. This observation can be demonstrated 
by examining the contrast evolutions of FBHs (or voids) with 
same characteristics but different widths (or diameters). 
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Thus, t 50 , is related to both the depth and diameter of the 
anomaly. During the time period from t 50i to t the sur- 
rounding area cools faster than the anomaly area. A larger 
diameter anomaly has a longer peak time. The peak time 
signifies the maximum in the lateral heat flow across the 5 
boundary of the trapped volume. The peak time is related to 
both the depth and the diameter. 

Beyond the peak time, the trapped volume continues to 
cool faster than the surrounding area. The lateral heat flow 
continues to reduce gradually from the peak time until the to 
contrast reaches a zero value. The lateral heat flow in the 
trapped volume dominates in the post peak time zone. The 
lateral heat flow is not only affected by the depth of the 
anomaly but also by the lateral transit time or the diameter of 
the anomaly. The t 5aR is also related to the anomaly width. 15 
The second inflection point may appear before or after the 
time t 5aR depending upon anomaly parameters. The second 
inflection point in the contrast evolution indicates transition 
from rapid cooling to slower cooling across the boundary of 
the trapped volume. 20 

To understand the overall direction of the heat flow from 
the measurement ROI, a term called the “effective heat flow” 
is introduced. The direction of the effective heat flow can be 
approximated as from the anomaly center point on the surface 
to the closest inner comer of the FBH. The tangent of this 25 
angle is given by y/2. The diffusivity in this direction can be 
called as an estimate of “the effective diffusivity” for the 
FBH. It is assumed that the effective diffusivity and the transit 
time from the center of the anomaly to the closest inner corner 
of the FBH have a controlling effect in influencing the con- 30 
trast evolution. The contrast prediction model uses y, which 
defines the direction of effective heat flow, to set-up the cali- 
bration. At any moment, the heat flow direction is different in 
different areas around the anomaly. For an anisotropic mate- 
rial, the diffusivity is different in different directions. For 35 
anisotropic materials, the diffusivity may also be different at 
different locations. Thus, in reality, the heat flows that influ- 
ence the contrast evolution may be subject to varying diffu- 
sivity and direction with time. The IR Contrast method 
assumes that the normal and lateral diffusivities do not differ 40 
significantly and can be approximated in a single value of the 
effective diffusivity, although some variation in the effective 
diffusivity with anomaly depth has been accommodated in 
the methods of the embodiments described herein. 

True delaminations are different from flat bottom holes. 45 
Unlike the circumstance of a FBH, heat flows behind the 
delamination. Moreover, the gap thickness between the two 
faces of the delamination may vary from the center towards 
the edge. The gap thickness may be largest in the center. From 
the center of the delamination to the edge, the gap may taper 50 
off to an intermittent contact/gap and finally to a zero gap or 
continuous material. The intermittent contact close to the 
edge provides a path for heat conduction. Both the heat con- 
duction through the edge and the heat conduction behind the 
delamination cause change in the direction of the effective 55 
heat flow, changing the contrast evolution and the IR Contrast 
prediction. In some cases, the delamination may have varying 
degree of intermittent contact resulting in a positive heat 
transmissivity. The degree of bridging affects the contrast 
evolution and the contrast prediction. 60 

In delaminations with the tapering gap, the measurement 
ROI is located close to the center of the anomaly where the 
gap thickness is highest. The contrast evolution is influenced 
by the effective heat flow direction that is primarily influ- 
enced by the thicker gap . The region bounded by the effective 65 
heat flow directions at the plane of the anomaly is called the 
primary region of measurement for the chosen measurement 
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ROI. Therefore, in this situation, the contrast analysis predicts 
the EFBH (or EUG) diameter that relates to the width of the 
delamination where the gap is relatively thicker. The bound- 
ary or width of the thicker gap can be assessed using the 
half-max technique, which will be described in further detail 
later herein. The thinner gap surrounding the thicker gap does 
not influence the contrast evolution appreciably, and the thin- 
ner gap region is considered to be outside the primary region 
of the measurement. A similar effect occurs when the delami- 
nation gap is intermittent such that there is bridging (or a thin 
gap) at several regions of a delamination. When the bridging 
occurs at multiple locations within a delamination, the effec- 
tive heat flow angle for the measurement ROI may be heavily 
influenced by the bridging. Localized high pixel intensity 
locations occur between the bridging locations. Any of the 
localized high pixel intensity locations can be chosen as the 
measurement ROI . Therefore, the primary region of measure- 
ment around the chosen measurement ROI may be much 
smaller than the overall width of the delamination. The half- 
max technique, depending upon the choice of measurement 
and reference ROI, may provide boundaries for the overall 
delamination or the primary region of measurement. 

The width of the anomaly indications, denoted by w, lm , can 
be directly measured from the IR data by choosing an appro- 
priate frame. The IR Contrast technique assumes that the 
delamination gap thickness is uniform. The pixel intensities 
along a scan line drawn across the linear indication or through 
the center of a circular indication are plotted against the pixel 
location. For a randomly oriented indication boundary, a scan 
line is established normal the boundary and the IR Contrast 
technique is used for localized edge detection. The maximum 
intensity pixel coordinate and its intensity (W max ) are noted. 
Similarly, the pixel coordinate and intensity of the baseline or 
reference level (W,.^ and W r ^) at left of the maximum 
intensity pixel and to the right of the maximum intensity pixel 
are noted. The average intensity of the maximum value and 
the base level value is used as the 50 percent or half-max 
intensity level (W ; , m ). Although, other percent boundaries 
(e.g. 30 percent boundary) are sometimes used to make the 
boundary estimates more or less conservative. The half-max 
algorithm can be applied to the normalized intensity data- 
cube. The normalized data-cube is obtained by subtracting 
the pre-flash frame intensity of each pixel from correspond- 
ing pixels of all post-flash frames (i.e. W"=W-W°). Equa- 
tions (47a) through (47e) provide various formulas to com- 
pute the half max intensity: 


W^, L + W max 

Whm,L = ^ 

(47a) 

Wr,f,R + 
W hm,R = 2 

(47b) 

u/ W re f,R + W re f,R+2W max 

4 

(47c) 

w ref + w max 

M'ftm = 2 .and 

(47d) 

Kef + Kax 
n _ rej max 

(47e) 


Two pixel coordinates, one at left 50 percent of intensity and 
the other at right 50 percent of intensity, are located. Pixels 
between the left 50 percent peak to the right 50 percent peak 
are counted and converted to length by multiplying the width 
of the image pixel to compute an approximate measure of the 
anomaly width. This technique is called the “width measure- 
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ment using the half-max technique.” The technique is applied 
to frames between t 50i and t k assuming that the surface 
texture effects occur before t 50£ . For a randomly oriented 
indication boundary, many scan lines are established along 
the entire boundary such that the lines are normal to the 5 
boundary. The scan lines are used to plot the pixel intensity 
along the scan line and locate the 50 percent intensity pixel as 
the location of the edge of the indication. A running contrast 
threshold (e.g. % C rthresho2d =\ percent) is used to compare 
maximum running contrast C rmax and screen the indication to 
data for noise. The expression for running contrast is shown in 
equation (470, 


C r, 


max 


w„ 


„ - Wr'f 
Kef 


X 1 00 > C r, threshold ■ 


(47f) 15 


The ratio of the normalized half-max intensity to the spatial 
noise shall typically be a value more than two for meaningful , 0 
boundary estimation, as expressed in equation (47g) 

(W hm "/W nol , e > 2). (47 g) 

Resolution of the half-max pixel can be improved by inter- 
polation within pixels. The half-max estimation is in error if 25 
the delamination thickness is not uniform. Therefore, the 
half-max width estimation should be complemented by the 
EFBH or EUG diameter estimations as needed to improve 
confidence in the estimation. FIG. 23(a) illustrates the half- 
max boundary extraction on selected flat bottom holes and 30 
slots. 

In order to compare the EFBH estimation with the 
observed width, a quantity called the “width ratio” is intro- 
duced. The width ratio (t|i, psi) is defined as the ratio of the 
simulation estimated EFBH (or EUG) diameter to the half- 35 
max calculated width as expressed in equations (48a) and 
(48b), 


, D efbh 
• 1 1 = — . 
w hm 

(48a) 

. Deug,%v 
<!>= — ■ — ■ 

(48b) 
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A correlation of the width ratio measurements is established 
on the calibration standard FBHs with gamma. The i|> ratio is 
observed to be close to one for round FBHs. If the anomaly is 
not like a gapping delamination such that heat is partially 
conducted or transmitted through the anomaly, then the 50 
EFBH diameter estimation is expected to be smaller. There- 
fore, t[' correlates to the degree of heat transmission through 
the anomaly. The heat transmission through the anomaly 
occurs due to heat conduction through the air present in the 
delamination gap, contact between the delamination faces, or 55 
the bridging of material in the delamination gap. It is hypoth- 
esized that a ratio close to one for round anomalies implies 
that the indication is similar to a flat bottom hole or a gapping 
delamination. The lower values of width ratio will have more 
inaccuracies in depth assessment. Therefore, another ratio 60 
called the amplitude ratio is introduced. The ratio can be 
calculated for linear indications of flat bottom slots or linear 
delaminations. Hie T) EFBH is expected to measure slightly 
more than the actual width due to effect of the length, but little 
effect is expected in the accuracy of w ; , m . Accordingly, com- 65 
paring i[' between a round FBH and a long slot with the same 
depth and width equal to the diameter of the FBH, a higher i|> 


should be expected for the slot. The effect of length contri- 
bution in the contrast evolution is minimal beyond a length of 
1 .5 times the width. 

In order to compare the EFBH diameter estimation with the 
measured half-max width, a quantity called the “amplitude 
ratio” is introduced. Similar to the width ratio, this ratio is 
intended to evaluate the transmissivity or the gap thickness of 
the anomaly. If a delamination has contact (0>T>- 1 ) between 
the faces, the heat will transmit through the delamination. A 
term related to the thermal mismatch factor (T) is multiplied 
to the exponential term inside the summation sign in equation 
(1). In the contrast evolution equation (23), it is assumed to 
have a value of minus one. A value of 0<T <- 1 results in lower 
peak contrast compared to the flat bottom hole of same size 
(T=- 1 ), but the overall shape and the peak time are advanced 
due to the change in direction of the effective heat flow (e.g., 
a decrease in angle 0). Here, the half-max width (w, !m ) is used 
as an input to the IR Contrast simulation. The depth and the 
amplitude constant inputs to the simulation are adjusted to 
match the simulation. The value of the diameter constant is 
obtained from the calibration. The amplitude constant is 
changed to provide a best match with the measured data. This 
amplitude constant is denoted by Ca. and is called the 
“anomaly evaluation amplitude constant.” The ratio of the 
evaluation amplitude constant and the calibrated amplitude 
constant is called the amplitude ratio and is denoted by (zeta), 
which is expressed in equation (49a), 



(49a) 


The amplitude ratio compares the peak contrast (amplitude) 
of the anomaly with the peak contrast of an equivalent flat 
bottom hole of same depth and diameter. Thus, a ratio close to 
one would imply that a round anomaly is similar to the FBH. 
A lower value would imply that there is material bridging, a 
thin gap, or contact within the anomaly. The lower values also 
are likely accompanied by mismatch in shapes of the mea- 
sured contrast and predicted FBH contrast evolutions. The 
circumstance of the material bridging provides adherence of 
the material above the anomaly, but a mere contact may not. 
The anomaly transmissivity measure or the contrast attenua- 
tion is defined in equation (49b), 

H-=K- (49b) 

The attenuation ranges from zero to one. Higher attenuation 
value indicates more heat transmission through the anomaly. 
A value of one indicates full transmission and a value of zero 
indicates no heat transmission, similar to the circumstance of 
no heat transmission in a FBH. The heat transmission through 
the anomaly is also due to partial contact or gap within the 
anomaly region. The heat may also transmit through the 
anomaly if the gap is too small, as thinner air gaps are more 
conductive than the larger gaps. Hie effect of anomaly trans- 
mission results in reducing the angle (0) of the effective heat 
flow. The effective heat flow is related to y and A c . Therefore, 
the calibration curves which are functions of A c can be devel- 
oped for anomalies with controlled transmission (e.g., known 
air gap thickness or effective gap conductivity and thickness). 
Thus, the EFBH analysis can be extended to assumed equiva- 
lent uniform gap (EUG) diameter and depth. 

A selected image area (or point) with a contiguous mosaic 
of fixed number of pixels is called the region of interest (ROI). 
The smallest ROI would contain just a single pixel. The next 
larger-sized ROI would comprise nine pixels in three rows 
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and three columns. Usually, the rows and columns of pixels 
are equal, with an odd number of pixel columns and rows 
providing symmetry about the center pixel. In some embodi- 
ments, the intensity of the ROI is the average of all pixel 
intensities contained in the ROT, and the pixel coordinate of 5 
the center pixel is assigned as the location of the ROI. Mea- 
surement accuracy of the half-max width is affected by the 
size of the ROI. Therefore, a term called the “pixel resolution 
ratio”, 2 is defined as expressed in equation (50a), 

^ (50a) 

R 

15 

In order to compare the ROI size with the width of the 
anomaly, another ratio parameter, X', is defined as expressed 
in equation (50b), 


where R is the size of the region of interest for the measure- 
ment pixel and w is the width or diameter of the anomaly. 25 
The size of a single pixel is also known as the instantaneous 
field of view (IFOV). A quantity called the slit response 
function (SRF) is determined experimentally for the IR cam- 
era to determine the resolution of the combination of the 
camera and lens. The slit response function provides a rela- 30 
tionship between the slit width (e.g. width of a warm object) 
at constant temperature with the measured relative intensity. 
Usually, the slit width should be a minimum of 3 to 5 pixels to 
register over 90 percent of maximum relative intensity. Thus, 
the width of an IR indication needs to be at least 4 times (X=4) 35 
the ROT size. An ROI with pixel size of 0.025 inch would be 
sufficient to measure the contrast on anomalies with width (or 
half-max width) larger than 0.1 inch. In other words, an 
indication with a half-max width of 0. 1 inch can be satisfac- 
torily measured with an ROT having a maximum dimension 40 
of 0.025 inch. Lower values of X cause averaging of the pixel 
intensity of the measurement ROI, resulting in lower contrast. 

If two ROIs of same size but different number of pixels are 
compared, the ROI with more pixels would have less noise in 
the intensity value. 45 

For curved surfaces, some defocusing is inevitable due to 
varying distances of the test object’s surface to the lens. The 
defocusing reduces the resolution by smearing the edges of 
the indication. In the above case, the ROI needs to be smaller 
than 0.025 inch to account for smearing due to the effect of 50 
defocusing. Although A values of larger than 4 may be satis- 
factory to measure contrast, it is not sufficient to measure the 
half-max width. The half-max width measurement resolution 
(ROI size) should be an order of magnitude smaller than the 
measurement (i.e., X>10). For X, values lower than ten are 55 
likely to introduce higher uncertainty in the half-max width 
measurement and the width ratio. Smaller pixel size, how- 
ever, provides better imaging of surface texture, but provides 
more variation in contrast data for detection of subsurface 
delaminations. Thus, a real small pixel size may not provide 60 
the most repeatable results. A pixel size that is one fourth of 
the feature size to be detected may provide a better compro- 
mise between the contrast repeatability and the resolution 
needed for the half-max width measurement. 

The peak contrast amplitude and peak contrast time are the 65 
two most important features that define the approximate 
shape of the contrast evolution curve. The simulation avail- 
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able from the ThermoCalc software program was used to 
study the correlation between the normalized contrast peak 
time and the peak amplitude for selected values of both the 
diameter and the depth of flat bottom holes. The simulation 
provided by ThermoCalc is based on finite difference mod- 
eling of the test object’s temperature, as shown in FIG. 20 . 
FIGS. 21 (a) and 21 (b) show similar contrast maps from the 
calibrated IR Contrast method described herein with respect 
to a test object made of RCC. The maps show similar trends 
but are different in magnitude due to the differences in mate- 
rial properties, assumptions, and accuracy of the models. The 
results of the simulation indicate that each combination of a 
FBH diameter and depth provides a unique combination of 
the corresponding peak time and peak amplitude for higher 
values of the contrast. Lines may be drawn through the con- 
stant diameter and the constant depth data points to form a 
grid. The constant diameter grid lines do not intersect for a 
peak contrast value above a certain threshold. The constant 
depth grid lines also do not intersect each other. The plot 
indicates that it is possible to assess the depth and diameter 
from the measured contrast evolution curve if the peak con- 
trast amplitude is above a certain threshold within the cali- 
bration space for the depth and diameter. The plot also sug- 
gests that the peak contrast is higher for higher values of D/d, 
e.g., y=20. The IR Contrast method is intended to capture the 
dependence of the peak contrast and the peak time on the 
diameter and depth using the calibration. The shape of the 
contrast evolution is usually taken from a thermally isotropic 
material with high y values. The IR Contrast prediction would 
be applicable to materials and anomalies with similar y val- 
ues. 

The IR Contrast method seeks to estimate the EFBH diam- 
eter and depth of the anomaly. As the peak contrast value 
decreases below a certain threshold, accuracy of establishing 
unique values of depth and diameter decreases. Accuracy in 
locating the peak contrast time and other times (e.g., t 50i and 
t 50A ) decreases significantly for lower peak contrast signals 
due to the correspondingly low signal to noise ratio. Deeper 
flaws provide long persistence time and lower peak ampli- 
tude, causing flattening of the contrast evolution curve at the 
peak point and reducing the accuracy of locating the peak 
time. Consequently, the accuracy of the depth measurement is 
reduced. Similarly, sensitivity of the peak contrast to the 
change in diameter decreases for deeper flaws. Thus, the 
accuracy of estimating y and the diameter also decreases for 
deeper flaws. Therefore, EFBH depth and diameter estima- 
tion is limited within a certain range of depth and y values. For 
deeper flaws, if an approximate value for one of the two 
parameters (e.g., depth or diameter) is known, the other may 
be assessed with better accuracy instead of estimating both 
the EFBH depth and diameter simultaneously. 

For example, consider a layered material with different 
thermal properties for each layer. The thermal response is 
affected by the variation in the layer thicknesses within a test 
object. The calibration standards should use the layer thick- 
ness covering the range of actual layer thicknesses. For a 
material with a coating layer, the thermal response is affected 
by the variation in thickness of the coating. An appropriate 
calibration file corresponding to the coating thickness should 
be used. The IR Contrast model uses three calibration con- 
stants (e.g., the amplitude constant, the diameter constant, 
and the depth constant) as the primary means to match the 
prediction to the measured data. Hie validation process 
should assess the degree of fit using any of the indices dis- 
cussed here. In some situations, it might be prudent to estab- 
lish the reference evolution from a reference anomaly (e.g. a 
FBH or flat bottom slot) in the desired material with nominal 
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properties and establish limits on prediction (y and at) in the 
neighboring region of the reference anomaly (e.g. y ref and 
M„y). The values of the evolution shape constants may also be 
changed in some embodiments. The model uses a single value 
for diffusivity for prediction of a single contrast evolution. 5 

Layered materials may be anisotropic in diffusivity. If the 
in-plane (lateral) and through thickness (normal) diffusivities 
are significantly different then the diffusivity experienced by 
the thermal conduction in the heat trapping volume, the sur- 
rounding area may be different at different post-flash times to 
and may cause the shape of the contrast evolution to be 
different compared to the contrast evolution from an isotropic 
material having the same value for y. Thus, the lateral heat 
flow in the trapped volume is affected by the lateral diffusiv- 
ity, while the heat flow in the surrounding area is dominated 15 
by the normal diffusivity. The contrast evolution shape is 
affected by interaction of these two heat flows. If lateral heat 
diffusivity is significantly different than normal diffusivity, 
then the contrast evolution curve beyond t peak is not likely to 
be useful for depth estimation, as it relates more to the differ- 20 
ence between the lateral diffusivity and the normal diffusiv- 
ity. The contrast evolution is dictated by heat diffusion in the 
reference and measurement ROIs. The primary factor affect- 
ing the contrast evolution near the peak time is heat diffusion 
in the direction dictated by y for the FBH. The diffusivity 25 
experienced by the heat flow not only depends upon the value 
of y, but also upon the layer thickness in the heat trapping 
volume and the in-plane and normal diffusivities in each 
layer. Thus, the diffusivity is dependent on the direction of the 
dominant heat flow. At different times in the contrast evolu- 30 
tion, the direction of the dominant heat flow changes, and 
therefore the effective diffusivity experienced by the heat 
flow can be different at different points of time during the 
contrast evolution for an anisotropic material. For an aniso- 
tropic material, correlation between the depth and the effec- 35 
tive diffusivity may be established to improve the contrast 
prediction. 

The depth estimation is likely to be more reliable based on 
the contrast evolution up to the peak time and particularly for 
the portion of curve between t 50i and t If the effective 40 
diffusivity does not differ significantly from the normal dif- 
fusivity, then the calibration profiles may compensate for the 
anisotropy. For a layered material with a given y, the effective 
diffusivity of the FBH may be dependent on the depth. 
Accordingly, in most embodiments the IR Contrast method 45 
needs a single diffusivity value and may be limited to a range 
of y values for an anisotropic material to limit the variation in 
the effective diffusivity. 

Comparison of the simulated evolution to the real evolution 
provides an estimation of the EFBH depth and diameter. 50 
Anomalies may or may not be circular in shape. Anomalies 
may be elongated in shape and have a portion of perimeter 
concave outward. In such situations, contrast is measured at a 
point of maximum localized intensity on a line drawn across 
the localized width of the indication. The intensity at this 55 
point is affected mostly by the closest boundary of the 
anomaly (e.g. width) and to a lesser extent by the continuation 
of the anomaly in the direction of length. The IR Contrast 
measurements on the elongated flat bottom holes indicate that 
as the aspect ratio (4>) increases beyond one, the length con- 60 
tributes to contrast amplitude, but the derivative of slope d 
C M /dl is negative, implying diminishing contribution of 
length to the contrast amplitude with increasing length. In 
practice, the contrast evolution is not affected by aspect ratio 
of 1 .5 and above for high gamma values. Consider a case with 65 
two anomalies at the same depth, one with a round shape, and 
the other with a square shape such that the diameter equals the 
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length of the side of the square anomaly. The square anomaly 
has 27 percent larger area compared to the round anomaly and 
consequently has a 13 percent larger equivalent FBH diam- 
eter. Thus, for elongated anomalies, the equivalent FBH 
diameter estimations are expected to be slightly larger than 
the physical width by 13 to 25 percent to account for the effect 
of the shape and length. The EFBH depth and diameter are 
expected to bear correlation to the actual width and depth of 
a delamination-like anomaly under certain limitations. FIG. 
23 (b) shows both the actual width and the equivalent diameter 
on holes numbered 3, 4, and 5 of an RCC test specimen as 
shown in FIG. 23(a). 

If cracks are present in the heat trapping volume, the ther- 
mal response or the contrast evolution is affected. Cracks may 
have partial contact between the faying surfaces. Depending 
on the degree of contact, gap thickness, size, location and 
number of cracks in the heat trapping volume and the sur- 
rounding material, the thermal response is affected. If the 
crack distribution (cracks per unit surface area) is uniform 
and the same as in the calibration standard, then the effect of 
cracks can be accomited for by calibration. Consider for 
example a special case of a linear anomaly with a single crack 
running lengthwise through the center, bisecting the heat 
trapping volume. The crack length is at least as long as the 
width. The measurement pixel is located in the center of the 
anomaly width and right next to the location of the crack. 
Therefore, as the heat conducts through the trapping volume, 
it reaches the crack boundary before it reaches the anomaly. 
The measurement pixel is located close to the crack, and the 
transit time from the measurement pixel to the crack is shorter 
than the transit time to the next closest boundary, which is at 
a depth of the delamination. The effect of the crack boundary 
is felt earlier by the measurement pixel. This effect results in 
the advancement of the peak contrast time or appearance of a 
separate peak. If the peak is advanced, it also may result in an 
increase in the peak contrast value. However, the heat contin- 
ues to diffuse to the anomaly and then around the edges of the 
anomaly. The persistence energy time is governed by the heat 
trapping volume and not by presence of a single crack. There- 
fore, as the contrast peak is advanced, the trailing side of the 
contrast evolution has lower amplitude so that the area under 
contrast evolution curve is more or less unchanged. The shape 
of the contrast evolution curve is also changed accordingly. 
The simulation based on the FBHs without cracks and match- 
ing EFBH simulation to the contrast evolutions from the 
crack containing anomalies result in a lower estimation of 
anomaly depth. The mismatch between the simulation and 
measured contrast evolutions at the peak contrast region may 
indicate the severity of the effect of the linear crack. 

The extraction of optimal contrast evolution from the IR 
data-cube helps to achieve full potential of the IR Contrast 
analysis. Typically, three shapes of delamination-like anoma- 
lies are observed: round, elliptical, and elongated in the form 
of a strip. Examining the case of a round anomaly, the mea- 
surement ROI is located close to the geometrical center as 
determined by the highest intensity pixel, which is a unique 
pixel and can be easily located through user interface or by a 
computer algorithm. The contrast evolution is influenced by 
the measurement and the reference ROIs. The reference ROI 
can be located at a chosen radial distance and at any clocking 
orientation from the measurement ROI. A discussion of locat- 
ing the measurement and reference ROIs follows: 

Radial Distance of the Reference ROI: 

The chosen radial distance is such that it locates the refer- 
ence ROI just outside the anomaly footprint where the pixel 
intensity has no influence from the anomaly. Moving the 
reference ROI closer to the measurement ROI, the contrast 
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amplitude starts reducing. If the reference ROI is located 
within the boundary of the anomaly footprint, it is likely that 
the contrast evolution has less amplitude than the maximum 
possible from the anomaly. Moving the reference ROI away 
from the anomaly, at a certain radial distance, the maximum 
in the contrast amplitude is reached. Moving further away, the 
reference ROI is less related to the measurement ROT and 
bias due to local (including lateral) heal flow effects at the 
reference ROI influencing the contrast evolution. Thus, there 
is an optimal radial distance range for the reference ROI 
where the contrast evolution is close to the ideal maximum 
possible from the anomaly and is mostly invariant with fur- 
ther increase in distance of the reference ROI. Typically, the 
location of the optimal reference ROI is determined by ana- 
lyzing a frame with time in a range t 5aR to t tail . The optimal 
reference ROI provides the shallowest possible depth estima- 
tion and the largest possible EFBH diameter estimation. 

Reference ROI Clocking on Round Anomalies: 

If there were no lateral heat currents (or temperature gra- 
dients), then the clocking of the reference ROI should not 
matter. However, in most circumstances, there are always 
some lateral heat currents due to the test object curvature, 
uneven flash heating, and change in material properties or 
thickness along the surface and proximity of the test object’s 
edge. The lateral heat current from the reference to the 
anomaly indicates that the reference area is relatively wanner 
and the contrast evolution will be lower than the tme contrast 
potential of the flaw. The reverse is true for the heat flow in the 
other direction. Assuming that the lateral heat flow (0 deg) has 
the same direction in the vicinity of the flaw where both the 
measurement ROT and the reference ROI data will be mea- 
sured, consider two reference ROIs, one in the upstream ( 1 80 
deg) of lateral heat flow and the other in the downstream (0 
deg) of the lateral heat flow with respect to the measurement 
ROI, both located at the same radial distance from the mea- 
surement ROI. The two reference ROIs have opposing effect 
on the contrast evolution, i.e., the upstream reference ROI 
causes the contrast evolution to lower, whereas the down- 
stream reference ROI causes it to be elevated. Thus, the act of 
averaging the two reference temperature profiles provides 
some compensation for the lateral heat flow. If the reference 
ROIs are clocked at right angles (+/-90 deg) to the heat flow, 
then the lateral heat flow has the least effect on the contrast 
evolution. Variation in the contrast evolution due to clocking 
can be minimized by choosing the reference ROI such that the 
contrast at the tail end of the contrast evolution (at t tail ) is 
nearly zero. Thus, a search of the optimum reference ROT 
providing zero tail end contrast is needed if a single reference 
ROI is used. Another method involves choosing many refer- 
ence ROIs clocked uniformly over the hill 360 degree rota- 
tion, each located at the same radial distance. The temperature 
profiles of these reference ROIs can be averaged to create an 
artificial reference ROI for contrast computation. In this case, 
it is not necessary to locate the optimal clocking, but only the 
optimal radial distance. The optimal clocking provides the 
lowest EFBH depth prediction and the EFBH diameter pre- 
diction somewhere in between the extremities of prediction 
from the un-optimized contrast evolutions. 

Variation in the Measurement ROI Location: The measure- 
ment ROI is chosen at a location where the peak amplitude (at 
t peak ) of the pixel intensity is a maximum across a scan line 
that is approximately normal to the anomaly boundary. The 
area of the measurement ROI is chosen to be small, i.e., it is 
the size of a single pixel to minimize the effect of the area 
averaging. Consider a case of a round FBH where the refer- 
ence ROI is located optimally and the location of the mea- 
surement ROI is changed from its optimal location. As the 
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measurement ROI moves away from the optimal location, the 
contrast peak amplitude reduces and the peak time increases 
resulting in the increase in the EFBH depth estimation. It is 
observed that the EFBH diameter estimation diminishes 
5 slightly. The optimal measurement ROI provides the shallow- 
est depth estimation. This observation can also be applied to 
a rectangular FBH. For instance, consider a rectangular FBH 
with the measurement ROI located at the center. The center 
location of the measurement ROI would be the optimum 
to location. Imagine a centerline that runs parallel to the length 
of the anomaly. Moving the measurement ROI along the 
centerline away from the optimum location, the contrast esti- 
mations are affected depending upon proximity of the refer- 
ence to the closest end of the anomaly. If the distance of the 
15 measurement ROI from the end is less than the width of the 
FBH, then a slight decrease in EFBH diameter is expected, 
similar to the case of a round FBH. Similarly, as the measure- 
ment ROI is moved closer to the end, an increase in EFBH 
depth is expected. The slight decrease in width estimation can 
20 be accounted as error in measurements . The increase in depth 
signifies a limit on the location of the measurement ROI. 
The Reference ROI Locations for Elongated Anomalies: 
For elongated and elliptical flaws, the uniform clocking 
over 360 degrees of rotation may not be applicable. Altema- 
25 tively, two reference ROIs that are located on either side of the 
anomaly may be possible, and averaging of the reference 
ROIs may provide some compensation to the adverse effect of 
lateral heat flow. The reference ROI optimization schemes 
work only under the assumption that nearly uniform lateral 
30 heat flow (surface temperature gradient) exists in the zone 
between the measurement ROI and the reference ROIs. Thus, 
although reference ROI optimization is possible, it may not 
provide full compensation, especially if radial distance is too 
large. Alternatively, one may locate the anomaly area in the 
35 center of the flash-hood where the flash intensity is most 
uniform and the test object surface can be most normal to the 
camera line of sight, resulting in the least amount of lateral 
heat currents at the center of hood. Measurements performed 
at the center of symmetry of lateral temperature gradients 
40 may result in better results, as the temperature gradient is 
lower in the center compared to any other area. 

Separate Reference ROI Evolution: 

In the previous discussion, it was assumed that the refer- 
ence ROI is located on the same part and same camera shot. 
45 However, there can be certain advantages of locating the 
reference ROI on a different part or different camera shot or 
both. The reference ROI evolution may be obtained from a 
separate part (e.g. NDE reference standard) in a “separate 
shot” taken on the reference standard. The reference standard 
50 may be located next to the part area being imaged so that part 
of the standard also appears in the “same shot”, and the 
reference ROI then can be located on the NDE reference 
standard. The reference standard can be chosen to have opti- 
mal properties including flatness. Similarly, if a separate shot 
55 is chosen for the reference ROI, the reference ROI may be 
located exactly at the location of measurement ROI, thus 
improving chances of illuminating the reference area with 
same intensity of illumination. However, these benefits may 
not be realized if the reference standard has different thermal 
60 properties, or the presence of a separate reference standard 
influences the heat flow within the part inspected. 

In some embodiments, better image resolution may be 
obtained using a small image pixel size. The measurement 
ROI is chosen at a location of the maximum intensity in the 
65 frame of maximum contrast (X peak ). It is located from the 
intensity versus the distance profile along a line drawn across 
the longer axis of the IR indication. Larger ROI size may 
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reduce the peak amplitude (peak clipping or averaging) due to 
averaging of intensity over many pixels. A X value of 4 may 
result in acceptable pixel intensity averaging. The reference 
ROI size is usually larger than the measurement ROI. The 
larger reference ROI provides better confidence in the tem- 
perature evolution of the reference area due to noise reduction 
and suppression of the texture effect. As the size of reference 
ROI is increased, its distance from the measurement ROI 
needs to be increased to prevent the reference ROI from 
running into the anomaly indication area. Typically, the ref- 
erence ROI shall be outside the anomaly footprint (as defined 
by the half-max measurement) at time larger than t 5aR . 
Increased distance between the measurement and the refer- 
ence ROI may also cause issues with averaging over larger 
variations in the lateral heat flow or surface temperature gra- 
dients. Thus, reference ROI size cannot be increased beyond 
a certain limit. Use of a zoom lens reduces the pixel size and 
can provide better X values. Stability of the camera with 
respect to the test object becomes more critical with a zoom 
lens due to the smaller pixel size. 

The anomaly indication is formed due to the heat diffusion. 
Immediately after the flash or at the beginning of the heat 
diffusion (t<t 50i ), the absorbed heat is concentrated at the 
surface. If the object surface has a surface texture in terms of 
surface roughness, then the incident heat flux has variation 
from peak to valley of the roughness topography. The peaks in 
the surface texture absorb higher flux compared to the valley 
areas, and the high points in the surface texture appear as hot 
spots with the size comparable to the physical size of each 
high point. The surface texture affects the lateral and normal 
heat conduction and is the other major factor causing hot 
spots in IR images immediately after flash time. A similar 
effect is observed when shallow surface cracks are present 
causing IR indication due to uneven cooling rates and the heat 
flow direction at either side of the crack. If a pixel size is larger 
than quarter of width (X<4) of these features, then the pixel 
intensity can capture only a portion of the actual spot inten- 
sity. For a pixel diameter of about 0.023 inch, the smallest IR 
indication half-max width can be about 0.092 inch in order to 
have a measurement above 90 percent of maximum possible 
intensity. Thus, smaller (area and amplitude) indications 
appearing in the earlier frames are not measured adequately. 
For example, if one were to use a zoom lens to change from 
one to five, a dramatic increase in the pixel intensity of texture 
indications would be expected where pixels are located at the 
hot spots, with a corresponding decrease in intensity where 
pixels are located at the cooler spots. Thus, if surface texture 
is present in the part, the earlier portion of the contrast evo- 
lution curve (t<t 50i ) may be more affected by the pixel size. 
With time, the heat diffuses and the texture indications fade 
away due to heat conduction. Smaller indications fade away 
before the larger ones. In most situations, indications due to 
subsurface delaminations and not due to surface texture are of 
interest. Delaminations larger than a certain threshold size are 
of interest. Therefore, the X anomaly f° r the threshold size 
should be greater than four. The apparent 'k a „ omah . measured 
at time t 50£ should be greater than four. The decay of the 
surface texture indications is dependent on many factors, such 
as the spatial frequency of the surface texture features and the 
depth of the surface features (including the surface crack 
depth). Thus, the texture indications are superimposed on the 
subsurface anomaly indication. FIG. 24 shows an early IR 
frame of flash thermography data for a graphite epoxy lami- 
nate. The image reveals a weave-like pattern due to surface 
texture with simple contrast evolutions illustrating the posi- 
tive and negative texture peaks. 
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Near surface resolution (depth) refers to the ability to 
resolve the challenge of shallow delaminations. Shallow 
anomalies have a shorter peak time and persistence. FIG. 25 
illustrates that a higher frame rate frequency provides better 
5 time resolution and definition of the contrast evolution. In the 
example of FIG. 25, the 60 Hz and lower frame rates do not 
capture the peak in the contrast evolution. In order to run 
automatic fit routines, one must capture the peak of the con- 
trast. The near surface indications may be affected by the 
to surface texture. Also, depending on extent of the texture 
effect, the near surface delaminations may not be discernible 
in the contrast evolution. A higher frame rate helps in improv- 
ing definition of the early portion of the contrast evolution. 
The near surface resolution is affected by the anomaly size (or 
15 gamma) and the pixel resolution, among other factors. 

The predicted contrast evolution curve is based on “cali- 
bration” of the numerical model to the real contrast evolution 
curves within certain limits on y and w. In most embodiments 
for contrast prediction, the material should be thermally iso- 
20 tropic, and the diffusivity and other thermal properties should 
be the same as the calibration reference standard. The cali- 
bration is usually done on a reference piece of the same 
material, and round flat bottom holes are used to simulate the 
anomalies. The reference evolution can be based on the ref- 
25 erence material itself. If the IR Contrast is calibrated on the 
flat bottom round holes, then the IR Contrast approach pre- 
dicts the diameter and depth of an equivalent flat bottom hole 
(EFBH) that provides a close match to the contrast measured 
on test hardware with a delamination-like anomaly. Compari- 
30 son of the predicted contrast evolution curve to the measured 
contrast evolution curve is done under the following assump- 
tions. The material properties between the calibration stan- 
dard and the test object inspected are the same or are corre- 
lated (e.g., thermal properties including coating thickness for 
35 a two-layer material). The test object’s thickness is similar to 
the calibration standard. The anomaly depth is much less than 
the test object’s thickness (e.g. co«l ). Values for omega that 
are close to one are not desirable, for in such situations the 
bottom of the test object would influence the contrast signa- 
40 ture. A precise upper limit on w can be validated through 
calibration. 

In the IR contrast method, it is assumed that the contrast 
evolution curve and measurements on the contrast evolution 
(the peak contrast, peak product time, persistence energy 
45 time, etc.) are only influenced by characteristics of the 
anomaly (e.g., anomaly or FBH depth, width, gap thickness 
and length) and not by any factors that are not related to the 
characteristics of the anomaly, such as proximity of the 
anomaly to the edge of the test object, proximity of the 
50 anomaly to the edge of the flash zone, a less than optimal 
choice of reference ROI location, lateral heat flow not related 
to the anomaly, and uneven flash heating. It is assumed that 
the contrast evolution curve is positive and has a shape similar 
to the contrast evolution curve for flat bottom holes (or other 
55 programmed anomalies) in the calibration reference stan- 
dard. Factors not related to the anomaly characteristics are 
assumed to be small and assumed not to distort the contrast 
evolution from its expected shape similar to an ideal shape 
(i.e., a positive skewed distribution with sharp rise and 
60 gradual decay). 

Natural anomalies are like delaminations and differ from 
FBHs. The contrast evolutions of the FBH and of delamina- 
tions of the same size are likely to differ to some extent due to 
additional mass below the delamination compared to the 
65 FBH. If the delamination-like controlled anomalies are avail- 
able for calibration, then it is advisable to establish the refer- 
ence evolution and the calibration curves to these anomalies 
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to improve contrast prediction. In general, the additional mass 
under the delamination should cool the surface above the 
anomaly faster, possibly implying the need to use slightly 
higher depth constants (e.g., a higher diffusivity ) in the simu- 
lation routine. 

The IR Contrast technique assumes a high emissivity value 
(such as by way of non-limiting example, a value >0.8). Tlius, 
the pixel intensity is assumed to be proportional to a positive 
power (3 to 10 approximately) of the absolute surface tem- 
perature. The polished inner sides of the flash-hood help to 
reflect the flash energy onto the test object, thereby maximiz- 
ing the incident energy on the test object. The lamps are hot 
after flash, and the test object reflects the heat radiation from 
the lamps in proportion to the test object’s reflectivity char- 
acteristic. The shiny inner sides of the flash-hood may act as 
four front-side mirrors and provide multiple mirror images of 
the flash lamps, resulting in a portion of pixel intensity that is 
not related to the test object’s surface temperature. If a test 
object has low emissivity, high emissivity paint can be 
applied to the test object to improve the emissivity to some 
extent. On the other hand, applying paint also introduces 
more variables, such as the thermal properties of the paint, 
introducing a potential new source of errors. 

The emissivity is dependent on wavelength of the incident 
radiation. The test object may have a different value of emis- 
sivity for visible and infrared wavelengths. Flash radiation 
comprises visible white light, but after the visible radiation 
stops the lamps radiate infrared radiation (part of afterglow) 
at much lower intensity. The slow decay after the peak inten- 
sity time in the flash is undesirable, as it has the effect of 
reducing contrast. A square flash pulse of a given duration 
provides the most contrast but is impractical. It is recom- 
mended to keep the flash duration constant (e.g. 3 ms) during 
calibration and inspection. In some embodiments, the test 
object should have high emissivity for both visible light and 
thermal radiation. The object surface angle with the camera’s 
line of sight affects thermal absorption, radiation, and the 
resulting contrast measurements. 

The flash-hood uses two long flash lamps with length along 
the inner roof and on either side of the camera lens. The lamps 
provide fairly uniform incident light radiation (<5 percent 
variation) at the hood opening. Thus, the heat illumination 
can be thought to have a direction parallel to “the center axis 
of the hood.” The camera lens is located at the center of the 
hood. In some embodiments, for a flat test object that is 
located at the hood opening with the camera “optics axis” 
(optic centerline) pointing to the center of hood opening, the 
line of sight from the camera lens center to the area at the 
center of hood is collinear with the optic axis and normal to 
the surface of the test object. In areas of the test object that are 
away from the center of the hood, the line of sight makes an 
angle called the line of sight angle with the normal to the 
surface of the test object. The heat emission in the direction of 
the line of sight is received by the camera. The line of sight 
angle increases as the distance of the area increases from the 
optic axis. 

For a diffuse radiator, the emissivity changes as a function 
of the cosine of the radiation angle measured from the line 
that is normal to the surface of the test object, and it is 
assumed to follow Lambert’s law. Thus, the measured inten- 
sity is of maximum value if the ROI is located on the optic 
axis. If both the measurement ROI and the reference ROI have 
about the same line of sight angle, then the effect of angle on 
the contrast is minimized. In at least some embodiments, 
centering the area where contrast measurements are extracted 
provides exemplary measurements. 


44 

In another embodiment, the test object is curved with a 
sloping step. In this embodiment, the direction of the surface 
normal is not collinear with the hood axis at many locations of 
the surface of the test object, particularly including the slop- 
5 ing surface. The angle between the surface normal and the 
hood axis can be considered to be the incident angle of the 
flash radiation. The intensity of the incident radiation is pro- 
portional to the cosine of the incident angle. The absorption of 
incident energy is also proportional to the cosine of the angle 
10 for a diffuse radiator. Tlius, surfaces with surface normals that 
are not collinear with the hood axis receive less flash radia- 
tion, and they absorb less of the received radiation. If the 
surface also makes a line of sight angle with the camera, then 
the emitted radiation is proportional to the cosine of the line of 
15 sight angle. Thus, the angle of the object surface affects the 
pixel intensity and measured temperature from the pixel 
intensity. The curvature of the test obj ect and the test obj ect’ s 
surface angle result in undesirable surface or lateral tempera- 
ture gradients that can be accommodated. 

20 The calibration profiles of the diameter and amplitude con- 

stants are about the same for thermally isotropic materials. 
This similarity implies that, once we calibrate the predicted 
evolution on a set of different FBHs with high y values, the 
calibration is applicable to other materials under certain con- 
25 ditions, such as the requirement of high emissivity and similar 
range of y values. The procedure involves calibrating the 
predicted contrast evolutions using the measured contrast on 
FBHs of a reference material with known diffusivity (a,. e/ ). 
The calibration involves plotting the diameter and amplitude 
30 constants (e.g. and C A ) versus the flaw size parameter 
(A c ) using the value of the depth constant equal to one (C d: 
ref= 1 ). To measure the diffusivity on an isotropic material of 
interest, FBHs with identical y values are machined in the 
sample. Measured contrast evolutions on the sample are 
obtained. Next, using the measured evolutions, the simulation 
fit with the closest actual diameter and depth as input is found. 
The diameter and amplitude constant from the calibration 
profiles of the reference material are used. The only variable 
in play while seeking the simulation fit is the depth constant, 
41 1 which is chosen by trial and error. The effective diffusivity of 
the material is determined using equation (51), 


C d ( 51 ) 

45 ~ & e,ref ^ ■> 

L d,ref 

For isotropic materials, the effective diffusivity is a measure 
of the material diffusivity. 

50 A flowchart for the process of measuring the effective 
thermal diffusivity using this method is provided in FIG. 26. 
At step 2605, flash thermography video data is recorded as 
digital frames for a known size flat bottom hole found in a 
material of interest. At step 2610, contrast evolution data on a 
55 flat bottom hole with known depth and diameter is extracted. 
At step 2615, the contrast evolution is computed, followed by 
the filtering of the contrast evolution at step 2620. At step 
2625, the peak contrast, the peak contrast time, the left fifty 
percent peak time, and the right fifty percent peak time are 
60 calculated. Next, at step 2630, a simulation contrast evolution 
match with known diameter and depth values and assumed 
depth constant values is performed. At step 2635, the delta 
square between the simulation contrast and the filtered con- 
trast is calculated. At decision step 2640, the calculated delta 
65 square is compared to the value from the previous iteration. If 
the change in value is less than a set specified value, then the 
effective diffusivity can be computed using equation (51) at 
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step 2645. If the change is greater than a set specified value, 
then the process continues by going to step 2650, where the 
depth constant value is changed, before looping back to step 
2630 to perform the step of matching the simulation contrast 
evolution to know values. 5 

The peak contrast and peak contrast time are two important 
parameters in the IR Contrast analysis. The peak contrast 
value is controlled by the size parameter Ac and it is indepen- 
dent of the diffusivity. For identical FBH in two materials 
with different diffusivities, the two peak contrast values of the to 
contrast evolutions are the same, but the peak contrast time 
for the material with the higher diffusivity is shorter com- 
pared to the other material. The effect of the increase in the 
diffusivity is similar to the effect of decreasing depth of the 
FBH, provided the value of y is kept about the same. The 15 
effect of doubling the depth is equivalent to changing diffu- 
sivity to about one fourth of the prior value. The effect of 
change in the diffusivity and depth is illustrated in FIG. 27. 
The change in the diffusivity of the test object due to mecha- 
nisms such as mass loss, off-gassing, or chemical reaction can 20 
be assessed by monitoring the peak time (or t 50i ). The change 
in peak time provides a measure to assess the change in 
diffusivity, provided that anomaly dimensions and transmis- 
sivity do not change. The lack of change in the peak contrast 
(including peak time) of an anomaly contrast evolution sug- 25 
gests that there is no change in the size of the anomaly and in 
transmissivity. 

The effect of change in the diameter of an anomaly as to the 
contrast is illustrated in FIG. 28. An increase in the diameter 
increases both the peak contrast and the peak contrast time. 30 
The effect of the anomaly’s transmissivity (or attenuation) on 
the contrast is illustrated in FIG. 29. Depending upon the 
anomaly transmissivity, the peak amplitude is affected. In one 
embodiment, the calibration curves for a selected attenuation 
can be generated with the help of an independent simulation, 35 
such as using the software program ThermoCalc. The 
approach first involves generation of the calibration curves 
corresponding to flat bottom holes. 

A flow diagram for generating calibration curves for flat 
bottom holes without attenuation is illustrated in FIG. 30. At 40 
step 3005, flash thermography video data is recorded as a 
series of digital frames for calibration using a set of flat 
bottom holes of known size and made of a material of interest. 
Next at step 3010, contrast evolution data on all flat bottom 
holes with known depth and diameter is extracted. At step 45 
3015, contrast evolution for selected contrast data in the hole 
of interest is computed and filtered. At step 3020, the peak 
contrast, the peak contrast time, the left fifty percent peak 
time, the right fifty percent peak time, and the flaw size 
parameter are calculated. At step 3025, simulation contrast 50 
evolution match with known diameter and depth values, 
depth constant value of one, and correct material diffusivity is 
performed as well as the selection of diameter and amplitude 
constant values. At step 3030, the delta square between the 
simulation contrast and the filtered contrast is calculated. At 55 
decision step 3035, the calculated delta square value is com- 
pared with the value from the previous iteration and if the 
change is less than a specified set value, then the diameter and 
amplitude constant for the flaw size parameter are finalized at 
step 3040. If not, the process continues to step 3052, wherein 60 
the diameter and/or the amplitude constant value(s) are 
changed, and then the process loops back to step 3025. After 
step 3040, the process continues to decision step 3045, 
wherein it is determined whether all calibration hole data of 
interest has been analyzed. If it has not, the process loops back 65 
to step 3015. If it has, then the process continues to step 3050, 
wherein the calibration data of diameter constant and ampli- 
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tude constant versus the flaw size parameter are plotted and a 
calibration data file is created. At step 3054, the calibration 
file is loaded into the IR contrast application. From here, a 
parallel set of steps is performed, including at step 3055 the 
creation of a contrast map and at step 3060 the calculation of 
the adjusted flaw size parameter as a function of gamma. 

Next, in another embodiment, a ThermoCalc simulation of 
the same size (diameter and depth) and different gap thick- 
nesses is undertaken. The gap thicknesses range from large 
(which are considered to be the same as the flat bottom hole) 
to thin gaps. The simulation provides the temperature contrast 
evolutions for each void of interest. Both the peak contrast 
and the peak time are noted for each contrast evolution. For 
each void, the peak time is divided by the peak time of the 
corresponding flat bottom hole (same diameter and same 
depth) to obtain the peak time ratio. For each void, the peak 
contrast is divided by the peak contrast of the corresponding 
flat bottom hole (same diameter and same depth) to obtain the 
peak contrast ratio or attenuation. A plot of the peak time ratio 
to the peak contrast ratio is generated as shown in FIG. 31(a). 
A table with values of the peak time ratio is generated at a 
selected attenuation (e.g. 50%) as shown in FIG. 31(6). 

FIG. 32 provides a flowchart of the process for generating 
attenuated calibration curves. The process begins at step 
3205, wherein flash thermography surface temperature 
response in thin gap delaminations and flat bottom holes are 
simulated using ThermoCalc or other simulation software. At 
step 3210, peak time ratio versus peak contrast ratio (or 
attenuation) plots are created for various gamma values for 
each gap thickness. At step 3215, the IR contrast simulation is 
opened and loaded with the calibration for flat bottom holes. 
The selected depth and diameter are used as inputs to obtain 
a simulated contrast, and the flat bottom hole peak contrast 
and peak time are noted. At step 3220, an attenuation level is 
chosen, for example p=50%, and an attenuated peak contrast 
is computed as u times the flat bottom hole peak contrast. The 
peak time ratio is then noted and the attenuated peak time is 
calculated as the peak time ratio times the flat bottom hole 
peak time. At step 3225, the diameter and amplitude constants 
are set (manually) to achieve correct values of attenuated peak 
contrast and peak time. The diameter and amplitude constants 
are noted as well as the corresponding flaw size parameter. At 
step 3230, a decision is made whether all relevant hole sizes 
have been analyzed covering the desired range for the flaw 
size parameter. If not, then the process returns to step 3215. If 
all relevant hole sizes have been analyzed, then the process 
continues to step 3235, wherein calibration curves of diam- 
eter and amplitude constant versus flaw size parameter for the 
chosen value of attenuation are generated and a calibration 
file is created. At step 3240, the attenuated calibration file for 
the attenuated contrast prediction is loaded in the IR contrast 
application. The process then continues by performing two 
additional steps. At step 3245, the measured filtered contrast 
evolution is loaded and either an automated or manual simu- 
lation match is performed. Also from step 3240, an attenuated 
contrast map is generated at step 3250. 

In some layered materials, the effective diffusivity may be 
depth dependent. If the diffusivity varies monotonically 
between a range of depths of interest, then calibration can 
account for the depth dependence of the diffusivity. The depth 
constant calibration is controlled by the empirical equation 
(52), 
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Cd,d ~ Cd,b + *S<fLoj 




(52) 


where C dd is the depth constant at depth d, C db is the depth 
constant for base depth (e.g. 2.9), is the slope for the depth 
constant (e.g., -0.9), d is the value for delamination, FBH, or 
uniform gap depth, d 4 is the base depth for a calibration FBH, 
and x=d/d 4 is the depth to base depth ratio or simply depth 
ratio. The amplitude constant may also be depth dependent. 
Equations (53) and (54) capture the relationship between the 
amplitude constant for base depth and variable depth, 


C hjd — C'n f d C A,b A and 


(53) 


15 


CAF,d = CaFM + 5V,^Lo: 
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where C^is the amplitude constant for a given depth d, C AF,d 
is the amplitude factor for a given depth d, C Ab ^ is the 
amplitude constant from the calibration curve for a base depth 
and given attenuation p, C AF.b is the amplitude factor at base 
depth, usually 1 , and S AFyd is the slope for the amplitude factor 2 s 
(e.g., -0.1). Moreover the calibration curves for given cali- 
bration can be customized for a chosen attenuation level using 
equations (55) and (56), 


~ c ,) t * and 

(55) 

Ca#l = Ca, cPa*i> 

(56) 


where p is the attenuation or transmissivity measure, C A(I is 
the diameter constant for attenuation p, C A0 (same as C D ) is 
the diameter constant for no attenuation, S 04i is the slope for 35 
the diameter constant, A c is the flaw size parameter, C ' A is the 
amplitude constant for attenuation p, C AJJ is the amplitude 
constant for no attenuation, is the slope for the amplitude 
constant and is chosen by best fit (for example r\ =0.25), which 
depends upon the value for p. FIG. 9(a) is a graph of the 40 
calibration curves for the diameter constant and amplitude 
constant. FIG. 9(b) shows an example of calibration curves 
for no attenuation (FBH) and with 30 percent and 50 percent 
attenuation. For thin gaps with 50 percent attenuation, the 
contrast map indicates much lower peak contrast and peak 45 
times. The contrast application can choose a FBH calibration 
or the corresponding percent attenuation calibration. For long 
indications, the EFBH or EUG width estimations are higher 
than the actual width due to the length to width aspect ratio of 
greater than one. If the data contains only long indications, 50 
one may set-up the calibration on long slot (aspect ratio >3) 
standards using width in place of the diameter. In this situa- 
tion, correction for the length effect will not be needed for 
most of the length of the indication except near the ends. 

Delaminations typically have their largest gap in the center, 55 
and the gap tapers off towards the edge or boundary of the 
delamination. There is a limit on the fractional range of 
attenuation (e.g., Ap ftm =0.5) that influences the contrast evo- 
lution, primarily depending upon the gap thickness and the 
gap thickness slope (dg/dw) profile in a tapering delamina- 60 
tion. The fractional range on the attenuation is defined in 
equation (57), 

(57) 

Here, the variable p. ml „ is the attenuation corresponding to the 65 
optimum location of the measurement ROI, which in most 
circumstances also corresponds to the location of the maxi- 
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mum gap thickness and therefore has a relatively minimum 
value. The term [i ed is the attenuation at a location close to 
the half-max boundary and defines the edge of the primary 
region of measurement. EFBH (or EUG) analysis assumes 
that the contrast evolution is within the prediction region as 
shown in the dashed lined box of FIG. 21(a). If the peak time 
and peak contrast coordinate pair lie outside the prediction 
region or on the boundaries of the prediction region, then the 
estimates will likely be inaccurate. A spatial pixel intensity 
profile along the major and the minor axis of the anomaly can 
be evaluated to determine if the anomaly falls within the 
contrast prediction region. Any shape other than a shape of 
sharp bell curve for the intensity profile along any of the major 
or minor axes would clearly indicate that the anomaly span 
exceeds the maximum prediction limit along that direction. 

An example of experimental data of width ratio (EFBH 
diameter/half-max width) versus the actual gap in the speci- 
mens made of RCC is provided in FIG. 33(a). The results 
show a correlation between the delamination gap and the 
width ratio. FIG. 33(A) shows an experimental correlation 
between a gap correlation parameter (peak product time/half- 
max width 2 ) and the gap thickness from the same flash ther- 
mography data used in FIG. 33(a). The results of FIG. 33(A) 
also show a correlation between the delamination gap and the 
gap correlation parameter. 

An example of simulation results using ThennoCalc are 
shown in FIG. 34. The simulation was done on uniform 
square voids. A gap correlation parameter is defined as the 
peak product time divided by the width ' 2 ’ 3 (i.e., the width to 
the power of 2.3). The simulation indicates a better correla- 
tion forthinner gaps. The simulation and experimental result s 
provide evidence of correlation between the gap correlation 
parameter and the gap thickness. 

The IR Contrast approach provides two methodologies to 
analyze the flash thermography video data for characteriza- 
tion of delamination-like anomalies. The first methodology 
involves measurements such as the peak contrast, the peak 
product time, the peak time to contrast ratio, the persistence 
energy time, and the peak time, which can be plotted against 
the location of the measurement. These measurements show 
some correlation to the anomaly’s depth (peak time to con- 
trast ratio), the width of the anomaly (half-max), and the gap 
of the anomaly (ratio of peak product time to half-max width) . 
These measurements are also applicable to measurement of 
surface cracks. An aspect of this approach is conversion of the 
flash thermography digital video data to the normalized con- 
trast datacube. The contrast values of each pixel and the 
extracted contrast evolution features (e.g. peak contrast, peak 
time) of each pixel are used to generate 2-D (frame) images. 

The second methodology involves matching the calibra- 
tion-based predicted contrast evolution to the measured con- 
trast evolution. The matching of the simulated evolution to the 
measured contrast evolution provides size and depth of the 
equivalent flat bottom hole or equivalent uniform gap with a 
known attenuation. The EFBH (or EUG) depth and diameter 
are shown to bear correlation to the width, the depth, and the 
gap thickness of the measured anomaly. Comparison of the 
EFBH (or EUG) evolutions with measured contrast evolu- 
tions, to provide evaluation of delaminations in terms of the 
depth (EFBH or EUG depth), the width (half-max), the gap 
thickness (width ratio of the EFBH or EUG diameter to half- 
max width) and the presence of surface cracks (EFBH diam- 
eter and depth estimations ), confirms correlation to the width, 
the depth, and the gap thickness of the measured anomaly. 

Another aspect of the first methodology is the half-max 
edge detection, which is considered to be the primary mea- 
surement of the width, and the IR Contrast EFBH (or EUG) 
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diameter prediction is considered to be the secondary. Any 
difference in the two width measurements is thought to be due 
to non-uniformity of the gap thickness, the value of the gap 
thickness, the heat transmissivity across the delamination, as 
well as the limitations of the two estimations. Both the half- 5 
max and the EFBH provide analysis of the area of delamina- 
tion within a certain fractional range of attenuation that is 
weighed heavily by the width containing a relatively thicker 
gap in the center (primary region of measurement) and is 
weighed hardly at all by the thin gap portion of the delami- to 
nation at the periphery. The width ratio provides correlation 
with the delamination gap thickness in the primary region of 
measurement. A flowchart that diagrams these methodolo- 
gies is given in FIG. 17. 

FIG. 35 is a block diagram illustrating an apparatus for 15 
characterizing anomalies in a material in accordance with 
embodiments described herein. Apparatus 3510 comprises a 
processor 3515 and a memory unit Processor 3515 is config- 
ured to perform computations and general control operations, 
and memory unit 3520 is configured to acquire data 3525, 20 
store a set of instructions to calibrate the IR camera system 
and process the data 3530, analyze the data 3535, and to 
evaluate an anomaly for size, depth, gap thickness, and mate- 
rial diffusivity 3540. 

In light of the principles and exemplary embodiments 25 
described and illustrated herein, it will be recognized that the 
exemplary embodiments can be modified in arrangement and 
detail without departing from such principles. Also, the fore- 
going discussion has focused on particular embodiments, but 
other configurations are contemplated. In particular, even 30 
though expressions such as “in one embodiment,” “in another 
embodiment,” or the like are used herein, these phrases are 
meant to reference the possibilities of embodiment generally, 
and are not intended to limit the embodiment described to one 
particular configuration. As used herein, these terms may 35 
reference the same or different embodiments that may be 
combinable into other embodiments. 

Similarly, although exemplary processes have been 
described with regard to particular operations performed in a 
particular sequence, numerous modifications could be 40 
applied to those processes to derive numerous alternative 
embodiments. For example, alternative embodiments may 
include processes that use fewer than all of the disclosed 
operations, processes that use additional operations, and pro- 
cesses in which the individual operations disclosed herein are 45 
combined, subdivided, rearranged, or otherwise altered. 

In view of the wide variety of useful permutations that may 
be readily derived from the exemplary embodiments 
described herein, this detailed description is intended to be 
illustrative only, and should not be taken as limiting the scope 50 
of the embodiments described herein. What is claimed as the 
invention, therefore, are all implementations that come within 
the scope of the following claims and all equivalents to such 
implementations . 

Thus, although only a few exemplary embodiments have 55 
been described in detail above, those skilled in the art will 
readily appreciate that many modifications are possible in the 
exemplary embodiments without materially departing from 
the novel teachings and advantages of the invention. Accord- 
ingly, all such modifications are intended to be included 60 
within the scope of the invention as defined in the following 
claims. In the claims, means-plus-function and step-plus- 
function clauses are intended to cover the structures or acts 
described herein as performing the recited function and not 
only structural equivalents, but also equivalent structures. 65 
Thus, while a nail and a screw may not be structural equiva- 
lents in that a nail employs a cylindrical surface to secure 
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wooden parts together, whereas a screw employs a helical 
surface, in the environment of fastening wooden parts, a nail 
and a screw may be equivalent structures. 

What is claimed is: 

1. A method for characterizing an anomaly in a material, 
the method comprising: 

extracting contrast data; 

measuring a contrast evolution; 

filtering the contrast evolution; 

measuring a peak amplitude of the contrast evolution; 

determining a diameter and a depth of the anomaly, 
wherein determining the diameter and the depth of the 
anomaly called direct depth estimation comprises: 
calculating a multiplier M and a flaw size parameter A c 
using an amplitude constant C A equal to one for a first 
iteration of determining the diameter and the depth of 
the anomaly and for subsequent iterations of deter- 
mining using an amplitude constant derived from an 
amplitude constant calibration curve; 
calculating an adjusted diameter to depth ratio y'; 
estimating the depth of the anomaly using a diameter 
constant C^, equal to one for the first iteration of deter- 
mining the diameter and the depth and for subsequent 
iterations of determining using a diameter constant 
derived from a diameter constant calibration curve; 
estimating the diameter of the anomaly; 
comparing the estimate of the depth of the anomaly after 
each iteration of estimating to the prior estimate of the 
depth of the anomaly to calculate a change in the 
estimate of the depth of the anomaly; and 

repeating the step of determining the diameter and the 
depth of the anomaly until the change in the estimate of 
the depth of the anomaly is less than a set value. 

2. The method of claim 1, wherein measuring a contrast 
evolution comprises calculating a normalized pixel intensity 
contrast, wherein calculating a normalized pixel intensity 
contrast comprises calculating a ratio of: 

a change in pixel intensity at a measurement region of 
interest minus a change in pixel intensity at a reference 
region of interest compared to the sum of the change in 
pixel intensity at the measurement region of interest and 
the change in pixel intensity at the reference region of 
interest. 

3. The method of claim 1, wherein measuring a contrast 
evolution comprises calculating a normalized temperature 
contrast, wherein calculating the normalized temperature 
contrast comprises calculating a ratio of: 

the difference between a change in temperature at a mea- 
surement region of interest and a change in temperature 
at a reference region of interest to the sum of the change 
in temperature at the measurement region of interest and 
the change in temperature at the reference region of 
interest. 

4. The method of claim 1, wherein filtering the contrast 
evolution comprises applying at least one of: an FFT filter, an 
infinite impulse response filter, a finite impulse response fil- 
ter, an exponential running average filter, an alpha beta filter, 
and a Kalman filter. 

5. The method of claim 1, wherein calculating a multiplier 
M comprises calculating a ratio of the peak amplitude of the 
contrast evolution that is measured to the peak amplitude of a 
reference contrast evolution. 
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6. The method of claim 5, wherein the flaw size parameter 
A c is calculated according to the equation of: 


A c = 


M 

e'ZC A C 5 ’ 


5 


wherein e' is a ratio of a normalized pixel intensity contrast to 
a normalized temperature contrast, £ is an amplitude ratio; 
and C 5 is a contrast evolution shape constant. 

7. The method of claim 1, wherein the amplitude constant 
derived from the amplitude constant calibration curve com- 
prises using the flaw size parameter A c calculated according 
to the equation of: 


Ac = 



wherein C 4 is a contrast evolution shape constant, and C 9 is a 
contrast evolution shape constant. 

8. The method of claim 1, wherein the adjusted diameter to 
depth ratio y' is calculated according to the equation of: 

25 


1 

r ' = (-(ln(l-A c ))C 4 )^, 


wherein C 4 is a contrast evolution shape constant, and C 9 is a 30 
contrast evolution shape constant. 

9. The method of claim 1, wherein estimating the depth (d) 
of the anomaly is calculated according to the equation of: 
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wherein d is the depth of the anomaly, t is a time; t' is a 
normalized time, a e is an effective thermal diffusivity, C, , C 3 , 
C 6 , C 7 , and C 8 are constants, and C d is a depth constant. 

10. The method of claim 1, wherein estimating the diam- 
eter (D) of the anomaly is calculated according to the equation 
of: 


yd 

: <V 


50 


wherein d is the depth of the anomaly. 

11. The method of claim 1, wherein the set value equals a 
user defined percentage. 

12. An apparatus for characterizing an anomaly in a mate- 55 
rial, the apparatus comprising: 

one or more processors; and 

one or more memory units coupled to the one or more 
processors, 

the apparatus being constructed and arranged to: 60 

extract contrast data; 
measure a contrast evolution; 
filter the contrast evolution; 

measure a peak amplitude of the contrast evolution; 
determine a diameter and a depth of the anomaly, 65 
wherein determining the diameter and the depth of the 
anomaly called the direct depth estimation comprises: 


calculating a multiplier M and a flaw size parameter 
A c using an amplitude constant C 4 equal to one for 
a first iteration of determining the diameter and the 
depth for subsequent iterations of determining 
using an amplitude constant derived from an ampli- 
tude constant calibration curve; 
calculating an adjusted diameter to depth ratio y'; 
estimating the depth of the anomaly using a diameter 
constant C^, equal to one for a first iteration of 
determining the diameter and the depth and for 
subsequent iterations of determining using a diam- 
eter constant derived from a diameter constant cali- 
bration curve; 

estimating the diameter of the anomaly; 
comparing the estimate of the depth of the anomaly 
after each iteration of estimating to the prior esti- 
mate of the depth of the anomaly to calculate a 
change in the estimate of the depth of the anomaly; 
and 

repeat the act of determining the diameter and the depth 
of the anomaly until the change in the estimate of the 
depth of the anomaly is less than a set value. 

13. The apparatus of claim 12, wherein the apparatus being 
constructed and arranged to measure a contrast evolution 
comprises means for calculating a normalized pixel intensity 
contrast evolution, wherein calculating the normalized pixel 
intensity contrast evolution comprises calculating a ratio of: 

the difference between a change in pixel intensity at a 
measurement region of interest and a change in pixel 
intensity at a reference region of interest to the sum of 
the change in pixel intensity at the measurement region 
of interest and the change in pixel intensity at the refer- 
ence region of interest. 

14. The apparatus of claim 12, wherein the apparatus being 
constructed and arranged to measure a contrast evolution 
comprises means for calculating a normalized temperature 
contrast evolution, wherein calculating the normalized tem- 
perature contrast evolution comprises calculating a ratio of; 

the difference between a change in temperature at a mea- 
surement region of interest and a change in temperature 
at a reference region of interest and to the sum of the 
change in temperature at the measurement region of 
interest and the change in temperature at the reference 
region of interest. 

15. The apparatus of claim 12, wherein the apparatus being 
constructed and arranged to filter the contrast evolution com- 
prises at least one of: an FFT filter, an infinite impulse 
response filter, a finite impulse response filter, an exponential 
numing average filter, an alpha beta filter, and a Kalman filter. 

16. The apparatus of claim 12, wherein the apparatus being 
constructed and arranged to calculate the multiplier M com- 
prises means for calculating a ratio of the peak amplitude of 
the contrast evolution that is measured to the peak amplitude 
of a reference contrast evolution. 

17. The apparatus of claim 16, wherein the flaw size param- 
eter A c is calculated according to the equation of: 


M 


wherein e' is a ratio of a normalized pixel intensity contrast to 
a normalized temperature contrast, £ is an amplitude ratio, 
and C 5 is a contrast evolution shape constant. 
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18. The apparatus of claim 12, wherein the amplitude con- 
stant derived from the amplitude constant calibration curve 
comprises using the flaw size parameter A c calculated accord- 
ing to the equation of: 

yc 9 
-<r 4 - 

wherein C 4 is a contrast evolution shape constant, and C 9 is a 
contrast evolution shape constant. 

19. The apparatus of claim 12, wherein the adjusted diam- 
eter to depth ratio y' is calculated according to the equation of: 

15 
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wherein C 4 is a contrast evolution shape constant, and C 9 is a 20 
contrast evolution shape constant. 

20. The apparatus of claim 12, wherein estimating the 
depth (d) of the anomaly is calculated according to the equa- 
tion of 
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wherein d is the depth of the anomaly, t is a time, t’ is a 
normalized time, a e is an effective thermal diffusivity, C t , C 3 , 

C 6 , C 7 , and C 8 are constants, and C d is a depth constant. 

21. The apparatus of claim 12, wherein estimating the 35 
diameter (D) of the anomaly is calculated according to the 
equation of: 


wherein d is the depth of the anomaly. 

22. The apparatus of claim 12, wherein the set value equals 45 
a user defined percentage. 

23. The apparatus of claim 12 further comprising: 

a flash lamp; 

an IR camera operatively connected to the one or more 
processors and the flash lamp; 50 

data acquisition electronics operatively connected to the IR 
camera; and 

a power supply/trigger unit operatively connected to the 
flash lamp, the power supply/trigger unit being opera- 
tively comiected to the one or more processors, the one 55 
or more memory units, and the data acquisition electron- 
ics, the one or more processors being configured with a 
set of instructions to trigger lighting of the flash lamp 
and to initiate the acquisition of the contrast data. 

24 . A computer program product stored on a non- transitory 60 
computer operable medium, the computer program product 
comprising software code being effective to characterize an 
anomaly in a material by comprising a set of instructions to: 

extract contrast data; 

measure a contrast evolution; 65 

filter the contrast evolution; 

measure a peak amplitude of the contrast evolution; 


determine a diameter and a depth of the anomaly, wherein 
determining the diameter and the depth of the anomaly 
called the direct depth estimation comprises: 
calculating a multiplier M and a flaw size parameter A^ 
using an amplitude constant C A equal to one for a first 
iteration of determining the diameter and the depth 
and for subsequent iterations of determining using an 
amplitude constant derived from an amplitude con- 
stant calibration curve; 

calculating an adjusted diameter to depth ratio y'; 

estimating the depth of the anomaly using a diameter 
constant C D equal to one for the first iteration of deter- 
mining the diameter and the depth of the anomaly and 
for subsequent iterations of determining using a diam- 
eter constant derived from a diameter constant cali- 
bration curve; 

estimating the diameter of the anomaly; 
comparing the estimate of the depth of the anomaly after 
each iteration of estimating to the prior estimate of the 
depth of the anomaly to calculate a change in the 
estimate of the depth of the anomaly; and 

repeat the step of determining the diameter and the depth of 
the anomaly until the change in the estimate of the depth 
of the anomaly is less than a set value. 

25. The computer program product of claim 24, wherein 
the software code being effective to measure a contrast evo- 
lution comprises a set of instructions to calculate a normal- 
ized pixel intensity contrast evolution, wherein the step of 
calculating the normalized pixel intensity contrast evolution 
comprises calculating a ratio of: 

the difference between a change in pixel intensity at a 
measurement region of interest and a change in pixel 
intensity at a reference region of interest to sum of the 
change in pixel intensity at the measurement region of 
interest and the change in pixel intensity at the reference 
region of interest. 

26. The computer program product of claim 24, wherein 
the software code being effective to measure a contrast evo- 
lution comprises a set of instructions to calculate a normal- 
ized temperature contrast evolution, wherein the step of cal- 
culating the normalized temperature contrast evolution 
comprises calculating a ratio of: 

the difference between a change in temperature at a mea- 
surement region of interest and a change in temperature 
at a reference region of interest to the sum of the change 
in temperature at the measurement region of interest and 
the change in temperature at the reference region of 
interest. 

27. The computer program product of claim 24, wherein 
the software code to filter the contrast evolution comprises a 
set of instructions for applying at least one of: an FFT filter, an 
infinite impulse response filter, a finite impulse response fil- 
ter, an exponential running average filter, an alpha beta filter, 
and a Kalman filter. 

28. The computer program product of claim 24, wherein 
the software code to calculate the multiplier M comprises a 
set of instructions for calculating a ratio of the peak amplitude 
of the contrast evolution that is measured to the peak ampli- 
tude of a reference contrast evolution. 

29. The computer program product of claim 28, wherein 
the flaw size parameter A c is calculated according to the 
equation of: 
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32. The computer program product of claim 24, wherein 
M the software code comprises a set of instructions for estimat- 
ed c A c 5 ’ ing the depth (d) of the anomaly according to the equation of 


wherein e 1 is a ratio of a normalized pixel intensity contrast to 
a normalized temperature contrast, ^ is an amplitude ratio, 
and C 5 is a contrast evolution shape constant. 

30. The computer program product of claim 24, wherein 
the software code being effective to use the amplitude con- 
stant derived from the amplitude constant calibration curve 
comprises a set of instructions for using the flaw size param- 
eter A c calculated according to the equation of: 



wherein C 4 is a contrast evolution shape constant, and C 9 is a 
contrast evolution shape constant. 

31. The computer program product of claim 24, wherein 
the software code being effective to calculate the adjusted 
diameter to depth ratio y’ comprises a set of instructions for 
calculating y' according to the equation of: 
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wherein C 4 is a contrast evolution shape constant, and C 9 is a 
contrast evolution shape constant. 
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wherein d is the depth of the anomaly, t is a time, t' is a 
normalized time, a e is an effective thermal diffusivity, C 1 , C 3 , 
C 6 , C 7 , and C 8 are constants, and C d is a depth constant. 

33. The computer program product of claim 24, wherein 
the software code comprises a set of instructions for estimat- 
ing the diameter (D) of the anomaly according to the equation 
of: 


wherein d is the depth of the anomaly. 

34. The computer program product of claim 24, wherein 
the set value equals a user defined percentage. 
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